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Col 1: IRAS name 

Col 2: Nuclear Separation in kpc unit 
Col 3: Tail length in kpc unit 
Col 4-' Morphological class 

Col 5: Spectral type: HII = H II region galaxy, L = LINER, S2 = Seyfert 2, SI = Seyfert 1 
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Parameters of Surface Brightness Profiles 
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OPTICAL AND NEAR-INFRARED IMAGING OF THE IRAS 1-JY SAMPLE OF ULTRALUMINOUS 

INFRARED GALAXIES. II. THE ANALYSIS 

S. Veilleux 1 - 2 , D.-C. Kim 3 - 4 , and D. B. Sanders 2,5 
ABSTRACT 

An R & K 1 atlas of the IRAS 1-Jy sample of 118 ultraluminous infrared galaxies (ULIGs) was presented in 
a companion paper (Kim, Veilleux, & Sanders 2002; Paper I). The present paper discusses the results from the 
analysis of these images supplemented with new spectroscopic data obtained at Keck. All but one object in the 
1-Jy sample show signs of a strong tidal interaction/merger. Multiple mergers involving more than two galaxies 
are seen in no more than 5 of the 118 (< 5%) systems. None of the 1-Jy sources is in the first-approach stage 
of the interaction, and most (56%) of them harbor a single disturbed nucleus and are therefore in the later stages 
of a merger. Seyfert galaxies (especially those of type 1), warm ULIGs (/25//6O > 0.2) and the more luminous 
systems (> 10 12 ' 5 Lq) all show a strong tendency to be advanced mergers with a single nucleus. 

The individual galaxies in the binary systems of the 1-Jy sample show a broad distribution in host magnitudes 
(luminosities) with a mean of -21.02 ± 0.76 mag. (0.85 ± L*) at R and -23.98 ± 1.25 mag. (0.90 ± ^ 4 
L*) at K' , and a R- or /if' -band luminosity ratio generally less than ~ 4. Single-nucleus ULIGs also show a broad 
distribution in host magnitudes (luminosities) with an average of -21.77 ± 0.92 mag. (1.69 ± H7 L*) at R and 
-25.03 ± 0.94 mag. (2.36 ± j'jg) at K'. These distributions overlap considerably with those of quasars. The same 
statement applies to R-K 1 colors in ULIG and quasar hosts. 

An analysis of the surface brightness profiles of the host galaxies in single-nucleus sources reveals that about 
73% of the R and K' surface brightness profiles are fit adequately by an elliptical-like /? 1 / 4 -law. These elliptical- 
like 1-Jy systems have luminosity and /?-band axial ratio distributions that are similar to those of normal (inactive) 
intermediate-luminosity ellipticals and follow with some scatter the same fi e - r e relation, giving credence to 
the idea that some of these objects may eventually become intermediate-luminosity elliptical galaxies if they 
get rid of their excess gas or transform this gas into stars. These elliptical-like hosts are most common among 
merger remnants with Seyfert 1 nuclei (83%), Seyfert 2 optical characteristics (69%) or mid-infrared (ISO) AGN 
signatures (80%). The mean half-light radius of these ULIGs is 4.80 ± 1.37 kpc at R and 3.48 ± 1.39 kpc at 
K', typical of intermediate-luminosity ellipticals. These values are in excellent agreement with recent quasar 
measurements obtained at H with HST, but are systematically lower than other HST measurements derived at R. 
The reason for this discrepancy between the two quasar datasets is not known. 

In general, the results from the present study are consistent with the merger-driven evolutionary sequence "cool 
ULIGs — > warm ULIGs — > quasars." However, many exceptions appear to exist to this simple picture (e.g., 46% 
of the 41 advanced mergers show no obvious signs of Seyfert activity). This underlines the importance of using 
a large homogeneous sample like the 1-Jy sample to draw statistically meaningful conclusions; the problems 
of small sample size and/or inhomogeneous selection criteria have plagued many studies of luminous infrared 
galaxies in the past. 

Subject headings: galaxies: active - galaxies: interactions - galaxies: Seyfert - galaxies: starburst - infrared: 
galaxies 



1. INTRODUCTION 

Local ultraluminous infrared galaxies (ULIGs; log [Lir/L q ] 
> 12; Ho = 75 km s" 1 Mpc" 1 and qo = 0) represent some of 
the best laboratories to study in detail the violent aftermaths of 
galaxy collisions and their possible connection with quasars and 
normal (inactive) elliptical galaxies. Recent deep surveys with 
the Infrared Space Observatory (ISO) and sub-mm ground- 
based facilities have revealed several distant (z « 0.5 - 4.0) 
infrared-luminous galaxies which appear to share many of the 
properties of local ULIGs (ISO : Kawara et al. 1998; Puget et 
al. 1999; Matsuhara et al. 2000; Efstathiou et al. 2000; Serjeant 
et al. 2001; Sanders et al. 2002, in prep.; SCUBA : Smail, Ivi- 



son, & Blain 1997; Hughes et al. 1998; Blain et al. 1999; Eales 
et al. 1999; Barger, Cowie, & Sanders 1999). This extragalactic 
population of high-z infrared bright galaxies appears to domi- 
nate the far-infrared extragalactic background and is probably a 
major contributor to the overall star formation and metal enrich- 
ment history of the universe (e.g., Smail et al. 1997; Hughes et 
al. 1998; Barger et al. 1998; Genzel & Cesarsky 2000; Rupke, 
Veilleux, & Sanders 2002; Blain et al. 2002). 

Two crucial questions need to be answered to properly ad- 
dress the issue of the origin and evolution of ULIGs: (1) What 
is the dominant energy source in ULIGs: Starbursts or active 
galactic nuclei (AGNs)? (2) Is the dominant energy source a 
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function of the interaction/merger phase in these systems? Con- 
siderable progress has been made in recent years in answering 
the first of these questions. Ground-based optical and near- 
infrared spectroscopic survey of the 1-Jy sample has shown 
that at least 25 - 30% of ULIGs show genuine signs of AGN 
activity (Kim, Veilleux, & Sanders 1998; Veilleux, Kim, & 
Sanders 1999a; Veilleux, Sanders, & Kim 1997, 1999b; see 
also Goldader et al. 1995; Murphy et al. 2001b). This frac- 
tion increases to 35 - 50% among the objects with log [L IR /L ] 

> 12.3. Comparisons of the dereddened emission-line lumi- 
nosities of the BLRs detected at optical or near-infrared wave- 
lengths in the ULIGs of the 1-Jy sample with those of opti- 
cal quasars indicate that the AGN/quasar in ULIGs is the main 
source of energy in at least 15 - 25% of all ULIGs in the 1- 
Jy sample. This fraction is closer to 30 - 50% among ULIGs 
with Lir > 10 12 - 3 Lq. These results are compatible with those 
from recent mid-infrared spectroscopic surveys carried out with 
ISO (e.g., Genzel et al. 1998; Lutz et al. 1998; Rigopoulou 
et al. 1999; Tran et al. 2001). Indeed, a detailed object-by- 
object comparison of the optical and mid-infrared classifica- 
tions shows an excellent agreement between the two classifica- 
tion schemes (Lutz, Veilleux, & Genzel 1999). These results 
suggest that strong nuclear activity, once triggered, quickly 
breaks the obscuring screen at least in certain directions, thus 
becoming detectable over a wide wavelength range. 

Much of the research effort now focusses on answering the 
second, more difficult question of a possible dependence of 
the energy source on the interaction/merger phase. A large 
dataset already exists in the literature on the morphology of 
luminous and ultraluminous infrared galaxies. Optical studies 
have shown that the fraction of strongly interacting/merger sys- 
tems increases with increasing infrared luminosities, reaching 

> 95% among the ultraluminous systems (e.g., Sanders et al. 
1988a; Melnick & Mirabel 1990; Murphy et al. 1996; Clements 
et al. 1996; although see Lawrence et al. 1989; Zou et al. 
1991; Leech et al. 1994). The improved angular resolution 
and lower optical depth in the near-infrared as compared to the 
optical provides a cleaner view of the nuclear stellar distribu- 
tion of ULIGs. The infrared morphologies of ULIGs often look 
significantly different from the optical images (e.g., Carico et 
al. 1990; Graham et al. 1990; Eales et al. 1990; Majewski 
et al. 1993). Over the past few years, high-resolution imag- 
ing with HST has contributed significantly to our knowledge of 
ULIGs (e.g., Surace et al. 1998; Zheng et al. 1999; Scoville 
et al. 2000; Borne et al. 2000; Cui et al. 2001; Farrah et al. 
2001; Colina et al. 2001; Bushouse et al. 2002). These stud- 
ies have had great success characterizing the central cores and 
massive stellar clusters in these objects, but the images often do 
not reach faint enough flux limits to fully characterize the host 
galaxies and associated tidal features. Adaptive optics imag- 
ing with ground-based telescopes offers a promising new way 
to acquire deep high-resolution of ULIGs, but this technique 
has so far been used for only a limited number of objects (e.g., 
Surace & Sanders 1999, 2000; Surace, Sanders, & Evans 2000, 
2001). 

In Kim, Veilleux, & Sanders (2002; Paper I), we presented 
an R & K' atlas of the IRAS 1-Jy sample of 1 18 ULIGs. This 
large and homogeneous dataset on the nearest (median redshift 
= 0.145) and brightest such objects in the universe is particu- 
larly well suited to address the issue of the origin and evolu- 
tion of ULIGs (see Kim 1995 and Kim & Sanders 1998 for a 
detailed description of the 1-Jy sample; note that Galactic ex- 



tinction is negligible for the 1-Jy sample since \b\ > 30° for 
all these sources). Results derived from the 1-Jy sample should 
serve as a good local baseline for studies on distant luminous 
infrared galaxies planned in the years to come with SIRTF, 
SOFIA, and other ground-based infrared and submm facilities. 
In the present paper, we analyze these images and combine the 
results of this analysis with those derived from our optical and 
near-infrared spectroscopic data to look for possible trends with 
optical and near-infrared morphological and spectrophotomet- 
ric parameters. When possible, the results from our analysis of 
the ground-based data on the 1-Jy sample are combined with 
published results obtained at other wavelengths and compared 
with the quickly growing set of high-quality data on optical and 
infrared-bright QSOs (e.g., McLure et al. 1999; Marquez et al. 
2000; McLeod & McLeod 2001; Canalizo & Stockton 2000a, 
2000b, 2001; Surace et al. 2001; Dunlop et al. 2002 and refer- 
ences therein; see Stockton 1999 for a review of the data before 
1999). Given the high frequency of mergers among ULIGs, the 
data on the 1-Jy sample also allow us to examine the disk-disk 
merger scenario for elliptical galaxy formation (e.g., Toomre 
1977; Schweizer 1982; Wright etal. 1990; Scoville et al. 1990; 
Stanford & Bushouse 1991; Kormendy & Sanders 1992; Doyon 
et al. 1994; Genzel et al. 2001). The results of our analysis are 
described in §§2 - 4. These results are then used in §5 to dis- 
cuss the origin and evolution of ULIGs and their possible evo- 
lutionary link with quasars and elliptical galaxies. A summary 
is presented in §6. Each object in the 1-Jy sample is described 
in detail in an Appendix. Preliminary results of this study were 
presented in Veilleux (2001). 

We adopt H = 75 km s" 1 Mpc" 1 and q = 0.0 throughout 
this paper and have converted the results from published papers 
to this cosmology to facilitate comparisons. We also compare 
our absolute magnitudes with the magnitudes corresponding to 
a L* galaxy in a Schechter function description of the local 
field galaxy luminosity function. The reader should be cautious 
when comparing luminosities from various papers as they may 
be based on different definitions of L* . The values of the abso- 
lute magnitudes of a L* galaxy in various bands, M*, are now 
better constrained than in the past thanks to large-scale surveys 
such as the 2dF Galaxy Redshift Survey (2dFGRS). Here we 
adopt M* R = -21 .2 mag. and M* K , = -24. 1 mag. The value of M* R 
was derived using a L* magnitude of M r = -20.9 mag. based on 
the results from the Las Campanas Redshift Survey (Lin et al. 
1996) and assuming a typical r-R = 0.3 mag. for a moderately 
old stellar population at z « 0.1 (Surace 2002, private commu- 
nication). A similar value for M R is derived using = -20.1 
mag. (e.g, Mobasher et al. 1993; Zucca et al. 1997; Ratcliffe et 
al. 1998; Folkes et al. 1999) and adopting B-R « 1.0 mag. for 
typical field galaxies at z ~ 0.1 (e.g. Fukugita et al. 1995), but 
keeping in mind that this color term is strongly dependent on 
morphological type and more specifically on the level of star 
formation activity in the galaxy. The K' magnitude is based 
directly on the results of Cole et al. (2001) from an infrared- 
selected subsample of the 2dFGRS. These results are also con- 
sistent to within ± <~ 0.2 mag. with the numbers quoted from 
optically selected samples (e.g., Mobasher, Sharpies, & Ellis 
1993; Loveday 2000) and from the K-band surveys by Glaze- 
brook et al. (1995) and Gardner et al. (1997), after taking into 
account the various k— and photometric corrections mentioned 
in Table 3 of Cole et al. (2001). 

2. GLOBAL AND NUCLEAR PROPERTIES 
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2.1. Absolute Magnitudes 

The integrated (nucleus + host) R and K' absolute magni- 
tudes and R-K' colors of the 1 -Jy ULIGs were discussed briefly 
in Paper I (§3.3). Figure 1 presents the global and nuclear (4- 
kpc diameter) color - magnitude diagrams of the 1-Jy sample. 
As mentioned in Paper I, the galaxies in our sample show a 
broad range of integrated R and K' luminosities ~ 0.4 - 120 
L* . There is no significant trend between global or nuclear R 
magnitudes and R-K' colors, although there is a slight ten- 
dency for the redder objects to have brighter K' magnitudes es- 
pecially when considering the nuclear quantities (Fig. 1). The 
more luminous 1-Jy objects at R and K' show a tendency to 
have warm IRAS 25 /im/60 /mi colors (/25//6O > 0.2) and to 
present Seyfert/AGN characteristics at optical and mid-infrared 
wavelengths (Fig. 2). There also appears to be a weak positive 
correlation between infrared luminosities and the global R or K' 
luminosities (see Figs. 2g and 2h). Note, however, that the 1-Jy 
sample is flux limited and is therefore subject to luminosity- 
dependent distance biases: systems with L 1R > 10 123 L Q are 
more distant on average (< z > = 0.183 ± 0.053) than ULIGs 
with L IR < 10 12 - 3 L Q (< z > = 0.132 ± 0.029). Consequently 
one should not be surprised to find trends between global lumi- 
nosities and infrared luminosities. No trends are seen between 
global R or K' luminosities and the IRAS 60 /im/100 /jm colors 
or the strengths (equivalent widths) of the stellar H/3 and Mg lb 
features. 

The ULIGs in the 1-Jy sample are on average less luminous 
at R than the quasars in the sample of Dunlop et al. (2002; 
see also Taylor et al. 1996; McLure et al. 1999), although 
there is considerable overlap between the two samples. Adopt- 
ing H - K 1 k, 0.8 for typical low-redshift quasars (e.g., PG 
QSOs; Surace 2002, private communication), we find that 1-Jy 
ULIGs are a better match in terms of total K' luminosities to the 
set of "low-luminosity" quasars of McLeod & Rieke (1994a; 
-23.2 < Mb < -22.1 mag. for our adopted Ho) than the "high- 
luminosity quasars" of McLeod & Rieke (1994b) and McLeod 
& McLeod (2001) with M B < -23.2 mag. However, we note 
in §4.2 that the situation is reversed when considering only the 
magnitudes of the host galaxies instead of the integrated mag- 
nitudes. 

2.2. Colors 

The (nuclei of the) 1-Jy ULIGs are (significantly) redder than 
normal galaxies [median R-K' and (R-K')4 = 3.25 and 4.27 
mag., respectively, versus 2.62 ± 0.34 mag. for normal galax- 
ies; see §3.3 in Paper I]. These colors are also redder than those 
of the PG QSOs (R-K' = 2.7 mag. on average if redshifted to 
Z « 0. 1 to match the average redshift of the 1-Jy ULIGs; Surace 
2002, private communication). The red ULIGs with integrated 
R-K' > 4 mag. or nuclear {R-K')^ > 5 mag. have a tendency 
to host Seyfert nuclei and have warm IRAS 25 /im/60 /im col- 
ors (Fig. 3). There is no significant correlation between the 
integrated or nuclear R-K' colors and the infrared luminosi- 
ties, IRAS 60 /mi/100 pm colors, and strengths of the stellar 
absorption features. Ten single-nucleus objects in the sample 
have (R-K')4 > 5 mag. and would therefore be considered 
extremely red objects (EROs; Elston, Rieke, & Rieke 1988) if 
it were not for the circumnuclear emission (which presumably 
would not be detected at slightly higher redshifts; neglecting 
fe-corrections). Of these ten objects, only two (Fill 19+3257 
and Fl 32 18+0552) have integrated colors which meet the ERO 
criterion; both of these objects are optically classified Seyfert 1 



galaxies. 

2.3. Compactness 

The compactness is defined as the ratio of the luminosity 
from the inner 4 kpc to the total luminosity. This quantity is 
physically meaningful only for single-nucleus systems. In Pa- 
per I, we mentioned that the mean compactness of our galaxies 
at K' (0.36±0.17) is significantly higher than at R (0.14±0.09). 
We find no significant trend between K' compactness and the 
infrared luminosity or IRAS feo / /100 color within the 1 -Jy sam- 
ple. However, Figure 4b shows that compact objects with 
(L K r)4/L K > > I host most of the optically identified Seyfert 
galaxies in our sample. Given the well-known correlations 
between warm IRAS 25 /im-to-60 /im color, the presence of 
an AGN, and small equivalent widths of Mg lb (especially in 
Seyfert Is; e.g., Veilleux et al. 1995, 1999a), it is not surpris- 
ing to find a tendency for these compact objects to also have 
/25//6O > 0.2 and EW(Mg lb) < 1 A. The trends with / 25 // 6 o 
and EW(Mg lb) are weaker than the trend with optical spectral 
type so it appears that the presence of a Seyfert nucleus is the 
key factor determining the value of the K' compactness. Similar 
trends are found when considering the R compactness. These 
results may reflect the fact that Seyfert nuclei contribute signif- 
icantly to the emission at R and K' or that the host galaxies of 
Seyfert nuclei are more compact than those of other galaxies. 
Our data on host galaxy properties appear to rule out this last 
possibility (see §4.5). 

3. MORPHOLOGICAL EVIDENCE FOR GALAXY 
INTERACTIONS 

3.1. Frequency of Multiple Mergers 

A careful inspection of the images in the 1-Jy atlas (Fig. 1 
in Paper I) combined with Keck spectroscopy of field galax- 
ies (described in detail in Appendix A) indicates that nearly 
100% of the objects in 1-Jy sample show signs of a strong 
tidal interaction/merger in the form of distorted or double nu- 
clei, tidal tails, bridges, and overlapping disks. The only possi- 
ble exception is F23233+2817 (see Appendix A). Nearly all of 
these objects are involved in the interaction/merger of two, and 
no more than two, galaxies. At most five systems out of 118 
(< 5%) are comprised of more than two interacting systems: 
F13443+0802, F14394+5332, F15001+1433, F17068+4027, 
and F21477+0502. In all five cases, triplets appear to be in- 
volved although the physical connection between the various 
components of the triplets have not all been confirmed. The 
fraction (number) of multiple (> 2) nuclei, double/pair merg- 
ers, and single-nucleus merger remnants in the 1-Jy sample is 
4% (5), 39% (46), and 56% (66, excluding F23233+2817), re- 
spectively. 

The general lack of multiple (> 2) systems in the 1-Jy sam- 
ple is not consistent with the relatively large fraction (~ 20%) 
of multiple systems claimed by Borne et al. (2000) and Cui et 
al. (2001) from an optical HST dataset. Part of the discrep- 
ancy may be in the sample selection. Very few objects in the 
HST sample are also part of the 1-Jy sample. It is possible, 
though probably not very likely, that the HST sample studied 
by both Borne et al. and Cui et al. just happens to include 
a larger fraction of multiple systems. Other sources of uncer- 
tainties will affect the measured fraction of multiple systems. 
The criteria for interaction differ slightly from one study to the 
next due to different observational parameters (e.g., high spatial 
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resolution will favor the detection of compact nuclei and stel- 
lar clusters; deep surface brightness limits will favor the detec- 
tion of faint tidal features) and the subjective nature of deciding 
if an interaction is indeed taking place without accurate spec- 
troscopy of all galaxies in the system. In our study, galaxies 
are said to be interacting only if they show convincing signs of 
large-scale tidal tails and/or present distorted outer isophotes. 
In many cases, Keck long-slit spectra of the galaxies in the sys- 
tems were obtained to confirm the interaction and to reject a 
few stellar sources or foreground/background galaxies. Some 
of the systems are accompanied by one or two galaxies at the 
same redshift as the IRAS sources but with no convincing signs 
of interaction. These galaxies are labeled "G" in Figure 1 of 
Paper I; they may be part of the same group of galaxies as the 
IRAS source but they do not appear to be interacting with it (see 
Appendix A). 

Another source of discrepancy is the definition of a galac- 
tic nucleus. Cui et al. define nuclei as having Mj < -17 mag. 
The HST study by Surace et al. (1998) of 12 warm ULIGs 
has found several genuine stellar superclusters with Mj < -17 
mag. These would be misclassified as nuclei according to the 
criterion of Cui et al. For comparison, all of the galaxies in 
the present paper have global Mr < -20. 1 mag. or Mj < -20.7 
mag. if we assume a typical R— I « 0.6 mag. (e.g., Fukugita et 
al. 1995). 

Finally, the waveband sampled by the observations is an 
important consideration when comparing the results from dif- 
ferent studies. The HST data in Borne et al. and Cui et 
al. were obtained in the /-band while our observations are 
taken at R and K' . The near-infrared data are particularly use- 
ful to reduce the effects of obscuration and thus to reveal the 
true stellar morphology of these galaxies. This is made obvi- 
ous when we compare our data with those of Cui et al. for 
the few objects that we have in common. Five of these ob- 
jects (F13539+2920, F14060+2919, F14202+2615, F22206- 
2715, F22491-1808) are part of the group of 17 multiple sys- 
tems of Cui et al. In all cases except F14202+2615, K' im- 
ages reveal the presence of only one or two galactic nuclei (in 
F14202+2615, the third galaxy shows no signs of interaction 
with the northern pair; see Appendix A). In all four of these 
cases, bright star-forming knots which were apparent in the R- 
band images disappear when observed at K', indicating the lack 
of an obvious old stellar population in these clusters. 

A recent HST //-band imaging study by Bushouse et al. 
(2002; see also Colina et al. 2001) of a subset of objects from 
the sample of Borne et al. (2002) and Cui et al. (2001) confirm 
our conclusions on the frequency of multiple systems. Out of 
27 ULIGs, Bushouse et al. find only one candidate for multi- 
ple merging (IRAS 18580+6527; see their Table 1), or - 4% of 
their sample. This fraction is virtually identical to that found 
from our sample. This excellent agreement emphasizes the 
importance of imaging galaxies at near-infrared wavelengths 
rather than in the optical when trying to determine the stel- 
lar distribution in these objects. We refer the reader to Ap- 
pendix A for a discussion of the four 1-Jy ULIGs (F11095- 
0238, F13469+5833, F20414-1651, and F22206-2715) which 
were also observed by Bushouse et al. (2002). 

In the rest of this section, we use the R and K' images of 
the 1-Jy sample to look for possible trends in the properties of 
ULIGs along the merger sequence. In the discussion below, 
three morphological parameters are used to trace the merger se- 
quence: the observed (projected) nuclear separation, the length 



of tidal tail(s), and the overall morphology of the system. 

3.2. Apparent Nuclear Separations 

Numerical simulations of merging disk systems indicate that 
the nuclear separation between the two galaxies in the system 
is not simply a monotonically decreasing function of time. The 
system often goes through a "pre-merger phase" (Mihos 1999 
calls it the "hanging out" phase) characterized by two distinct 
galaxies with well-developed tidal tails and bridges as a result 
of the first close encounter (e.g., Barnes & Hernquist 1992; Mi- 
hos & Hernquist 1996; Gerritsen 1997; Bekki et al. 1999). 
This result complicates the interpretation of the measured (pro- 
jected) nuclear separations (NS) in the 1-Jy sample of galaxies 
since, without kinematic information, we cannot tell if the two 
galaxies are pre- or post-apgalacticon (farthest approach). The 
measurements are listed in Table 1 and shown in Figure 5. An 
excellent agreement is found between our measurements and 
those of Surace (1998; adjusted to qo = 0). The distribution of 
nuclear separations is highly peaked at small values but also 
presents a significant tail at high values. The median separa- 
tion of the sample is below the seeing limit (< 1" ~ 3 kpc). 
More than 70% (60%) of the objects in the 1-Jy sample have 
nuclear separations of less than 5.0 (2.5) kpc, but 16% (7%) of 
the objects have separations in excess of 10 (20) kpc. These 
widely separated pairs may be difficult to explain in the stan- 
dard ULIGs scenario, where the main activity is triggered in the 
late phase of the galaxy merger (e.g., Mihos & Hernquist 1996), 
unless the ULIG phase was triggered by an earlier merger event 
which is not evident at the present epoch. This isssue is dis- 
cussed in more detail in §5.1. 

We have looked for possible correlations between the mea- 
sured NS and the optical spectral type, infrared luminosity, 
IRAS 25 pm/60 p,m and 60 /im/100 p,m colors within the 1-Jy 
sample. We find a significant tendency for the Seyfert 1 galax- 
ies and, to a lesser extent, the Seyfert 2 galaxies to have small 
nuclear separations. This is shown in Figure 6a. 100% (68%) 
of the Seyfert Is (2s) have NS < 2.5 (5.0) kpc. Similar trends 
are observed with IRAS fis/ feo colors (Fig. 6b): 79% of the 
14 warm objects with f 2 s/ fa) > 0.2 have nuclear separations 
of 2.5 kpc or less (all of them have a nuclear separation of less 
than 10 kpc). 

The trend between NS and AGN activity is confirmed when 
we use the ISO spectral classification of Lutz et al. (1999). 
Twenty-five objects are in common with our sample. Of the 
10 objects classified as AGN-dominated according to ISO, all 
of them have NS < 2.5 kpc. In contrast, starburst-dominated 
galaxies show a broad range of nuclear separation. 

No obvious trends are observed between the measured nu- 
clear separations and the equivalent widths of the H/3 and Mg 
lb stellar absorption features, two traditional indicators of the 
starburst age (Veilleux et al. 1995, 1999a). 

As shown in Figure 7, single-nucleus sources and close pairs 
are more common in the highest infrared luminosity bin (83% 
of the sources with log [Li R /L Q ] > 12.5 have NS < 2.5 kpc ver- 
sus 58% at log [Lir/L ] < 12.24). This trend appears to extend 
to lower infrared luminosities: In a recent study of 56 galax- 
ies with log [Lir/L ] = 11.10 - 11.99, Ishida (2002) finds that 
only about a third of the objects with log [Lir/L q ] > 11.5 are 
advanced mergers with a single nucleus. 

3.3. Lengths of Tidal Tails 
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The strength and length of tidal tails in galaxy mergers are 
strong functions of the encounter geometry and the merger 
phase. A strong resonance between the orbital and rotational 
motions of stars and gas produces strong tidal tails in prograde 
disks (e.g., Toomre & Toomre 1972; Mihos & Hernquist 1996). 
The lack of resonance in retrograde disks inhibits the formation 
of extended tidal tails, but still results in the formation of tidally 
induced spiral arms, central bars, and diffuse tidal debris (e.g., 
Mihos & Hernquist 1996). Systems in the pre-merger phase, 
where the two galaxies are still distinct but have gone through 
the first encounter, are expected to have well-developed tidal 
features. As the final merger takes place, the tidal tails gradu- 
ally disappear as the material in the tails is accreted back onto 
the remnant or escapes the system altogether (e.g., Mihos, Du- 
binski, & Hernquist 1998; Hibbard & van Gorkom 1996). 

In an attempt to further constrain the phase of the interaction 
in 1-Jy ULIGs, we have measured the total projected lengths of 
the tidal tails from our /?-band images (the K' images are less 
sensitive to these faint features, so they were not used for this 
analysis). The measurements were done along each tail down 
to a constant R surface brightness of 24 magnitudes per square 
arcsecond. The results are listed in Table 1 and shown in Figure 
8. Slightly more than half of the objects in the 1-Jy sample have 
tail lengths of 10 kpc or less, indicating an advanced merger. 
There is no obvious correlation between tail lengths and optical 
spectral types, infrared luminosities, and IRAS 25 ^m/60 ^m 
and 60 ^m/100 ^m colors. Note, however, that 70% of the ob- 
jects with /25 /f(,o > 0.2 have tail lengths of 10 kpc or less. This 
fraction is 60%, 36%, 52%, and 66% when considering Seyfert 
Is, Seyfert 2s, LINERs, and H II galaxies, respectively. 

3.4. Interaction Classification 

Using a classification scheme first proposed by Surace ( 1 998) 
and based on the results of published numerical simulations 
(e.g., Barnes & Hernquist 1992; Mihos & Hernquist 1996; Ger- 
ritsen 1997), we have classified each object in the 1-Jy sample 
(except for the five triplet candidates and F23233+2817, which 
is an isolated spiral galaxy) according to their morphology: 

I. First Approach. This phase refers to the earliest stage 
of the interaction, prior to the first close passage of the 
galaxies, when the galaxy disks remain relatively unper- 
turbed and separate and hence do not show evidence for 
tidal tails or bridges. 

II. First Contact. At this stage, the disks overlap but strong 
bars and tidal tails have not yet formed. These are the 
primary morphological differences between this phase 
and the early merger phase (Phase IV below). 

III. Pre-Merger. This stage is characterized by two identifi- 
able nuclei with well-developed tidal tails and bridges. 
While our imaging observations cannot prove that the 
double-nucleus systems will necessarily merge, the pres- 
ence of ultraluminous activity combined with the preva- 
lence of mergers in ultraluminous galaxies, as well as the 
presence of tidal tails and other structure interconnecting 
the nuclei in every case, would seem to strongly indicate 
eventual merger. Within this category the systems were 
further divided into two subclasses based on their pro- 
jected separations: 

Ilia. Wide Binary. Systems with apparent separation > 
10 kpc. 



Illb. Close Binary. Systems with apparent separation < 
10 kpc. 

Note that these subclasses are an approximation only to 
the merger phase since without detailed kinematic infor- 
mation for many of these systems, we cannot tell if the 
two galaxies are pre- or post-apgalacticon (farthest ap- 
proach). 

IV. Merger. This stage occurs after the nuclei have appar- 
ently coalesced. These systems have prominent tidal fea- 
tures, but only one nucleus can be detected at optical and 
near-infrared wavelengths. Additionally, the sole galaxy 
core is often noticeably extended and tends to be cut by 
many dust lanes. Following Surace (1998), we further 
subdivide this class into two categories: 
IVa. Diffuse Merger. Systems with (L K ')4/L K > < |. 
IVb. Compact Merger. Systems with (L K ')4/L K ' > |. 
Objects in class IVa consists of systems that have diffuse, 
extended central regions typically composed of several 
smaller extended emission region which appear to be cut 
by dust lanes. Objects in class IVb are dominated, partic- 
ularly at long wavelengths (K'), by a single point source. 
As pointed out by Surace (1998), there is no a priori rea- 
son to believe that the stage IVb systems are any older 
than the stage IVa systems, other than the fact that stage 
IVa resembles stage III, and stage IVb resembles stage 
V, and that stage V is older than stage III based on the 
appearance of tidal structure. 

V. Old Merger. These are systems which do not show 
any direct unmistakable signs of tidal tails, yet have 
disturbed central morphologies similar to those of the 
merger stage IV systems, particularly those with knots 
of star formation. At this stage, the surface brightness of 
the tidal features in the system has fallen below our de- 
tection limit, leaving nothing but the high surface bright- 
ness relaxed merger remnant core visible. 

The results of the classification are given in Table 1. The 
distribution is shown in Figure 9. None of the 1-Jy sources 
appears to be in the early stage of the interaction (I and II). Sys- 
tems comprised of two distinct objects all show prominent tidal 
features. Most (56%) of the ULIGs harbor a single nucleus and 
are therefore in the later stages of the merger. This is a very sig- 
nificant result: the ULIG phase is triggered only after the first 
close encounter between the two galaxies and generally when 
the two galaxies have coalesced into a single object. The frac- 
tion of singles to doubles in the 1-Jy sample varies with infrared 
luminosity, increasing steeply for the most extreme ULIGs with 
log [L IR /L Q ] > 12.5 (Fig. 10a). This fraction seems to decrease 
when considering lower luminosity systems with log [Lir/L©] 
= 1 1 .00 - 1 1 .99 (Ishida 2002). This class of objects is generally 
comprised of two distinct galaxies often widely separated and 
sometimes showing no obvious tidal features. Many of these lu- 
minous infrared galaxies (i.e. LIGs: log[Li R /L Q ] =11- 1 1.99) 
therefore fall into interaction classes I - III. 

Figure 10b shows that the interaction class is a strong func- 
tion of the optical spectral type. We find that all ten Seyfert Is 
are advanced mergers (class IVb or V) and that half (50%) of 
the Seyfert 2s are early or advanced mergers (class IVa, IVb, or 
V). This is consistent with the results of our analysis of their 
nuclear separations (§3.2). On the other hand, LINERs and 
H II region-like galaxies show no preference in interaction class 
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(other than to avoid classes I and II). A strong correlation is also 
found when we consider the warm objects with /25//6O > 0.2 
(Fig. 10c): Nearly 80% of these systems (11 objects out of 14) 
are advanced mergers (class IVb or V). 

We find no obvious trend between the interaction class and 
IRAS 60 i«n/100 /zm color. The strengths of the H/3 or Mg lb 
stellar features measured in the nuclei of the 1-Jy objects (either 
considering all the objects in the sample or only in systems op- 
tically identified as starburst galaxies) show no correlation with 
interaction class. This negative result indicates that traditional 
nuclear indicators of the starburst age are not particularly good 
indicators of the merger phase or the epoch of the merger event. 
Detailed studies of star clusters in merger-induced starbursts 
confirm the existence of a broad range of stellar age which re- 
flects the complex star formation history of these systems (e.g., 
Surace et al. 1998; Whitmore et al. 1999; Gallagher et al. 2001; 
Zhang, Fall, & Whitmore 2001, and references therein). 

4. PROPERTIES OF THE HOST GALAXIES 

4. 1 . Surface Brightness Profiles 

In an attempt to characterize the brightness distribution of 
the host galaxies in the ULIGs of the 1-Jy sample, we have per- 
formed the following isophotal analysis. First, a central point 
source was removed from the centers of the objects with obvi- 
ous central excess emission. For this procedure, a stellar source 
within the field of view of the detector was shifted and scaled 
appropriately in intensity to produce a smooth and monotonic 
profile in the central regions of the galaxies. Ellipses were then 
fitted to the residuals using the standard Fourier expansion (e.g., 
Binney & Merrifield 1998) 

4 

5((j)) = S + ^[A„ cos(n0) + B„ sin(n</>)], (1) 

n=l 

where 8(cf)) = Ri(<p) - R e (4>) is the distance at position angle <f) 
between the radii of the corresponding points on the ellipse and 
on the isophote, and A„ and B n are the coefficients of the Fourier 
expansion. The typical error for the ellipse fitting is estimated 
to be less than 0.1 mag. arcsec" 2 . Unless otherwise noted, the 
error bars on the various quantities derived from these fits are 
equivalent to 1-sigma values and were directly derived from 
the IRAF routines (see the IRAF help page and in particular 
Jedrezejewski 1987 and Busko 1996 for a more detailed dis- 
cussion of the error analysis). 

The R and K' surface brightness profiles derived from the el- 
lipse fitting (not corrected for cosmological dimming) are pre- 
sented in Figure 1 1 using two different scales on the horizontal 
axis: a logarithmic scale on the left and a linear 1 / 4 scale on the 
right. An exponential disk and an elliptical-like object with a 
de Vaucouleurs profile would produce straight lines in the left 
and right panels, respectively. A vertical arrow in each panel 
represents the size of the seeing disk. Least-square fits to the 
profiles were carried out to determine whether an exponential 
or de Vaucouleurs profile best fits the observed profile at R and 
K' . To avoid complications due to the seeing, the least-squares 
fit was applied to the data points only at radii larger than twice 
the seeing disk. The value of the reduced chi-square is indi- 
cated in the upper right hand corner of each panel in Figure 
1 1 . The results from the fits are summarized in Table 2. In the 
calculation of the reduced chi-square, the number of degrees of 
freedom is the number of data points in the observed radial pro- 
file minus the number of free parameters in the fit (= 2 for both 



the exponential and de Vaucouleurs profiles). A reduced chi- 
square of order unity indicates a good fit. If the chi-square is 
much larger than unity, then the fit is poor. If the reduced chi- 
square is much less than unity, then there is not enough data 
points to carry out the fit and the parameters derived from the 
fit are not reliable. 

Figure 12a shows the fraction of single-nucleus objects (i.e. 
merger remnants with interaction class IV or V) in the 1-Jy sam- 
ple with pure disk-like (D: 2%), pure elliptical-like (E: 35%), 
elliptical/disk-like (E/D: 38%), and ambiguous (A: 26%) radial 
profiles. "E/D" profiles are defined as being equally well fit 
by an exponential or de Vaucouleurs profile, while neither fits 
galaxies with "ambiguous" profiles. Both E/D and A profiles 
are quite common among ULIGs. The large number of E/D 
galaxies in our sample probably reflects the limitations of our 
data: images with a larger dynamical range in surface bright- 
ness (i.e. higher spatial resolution to probe the nuclear region 
and fainter surface brightness limits to map the outer portions 
of the galaxies) would help rule out one type of profile over 
the other. The large fraction of ambiguous profiles probably re- 
flects the residual effects of the galaxy interaction on the mor- 
phology of the merger remnants. If we combine the E and E/D 
galaxies, we find that 73% of the single-nucleus objects are fit 
adequately by a de Vaucouleurs profile. 

The fraction of ambiguous profiles among double-nucleus 
systems (Fig. 12b) is significantly larger than among single- 
nucleus systems. Several binaries show obvious tidal tails and 
bridges which distort the radial profile. Among the other bi- 
naries, we find that the fraction of elliptical-like galaxies is 
slightly larger than disk-like hosts, perhaps an indication that 
the progenitor galaxies generally have strong spheroidal com- 
ponents (e.g., bulges). 

Strong correlations are found between the type of host galaxy 
profile and the nuclear and global properties of single-nucleus 
systems. In the following analysis, we take the conservative ap- 
proach to double-count each E/D object as an E object and a D 
object (it is conservative because it has the effect of smoothing 
out possibly stronger trends with host galaxy profiles). Figure 
13a shows that the hosts of most merger remnants with Seyfert 
1 (83%) or Seyfert 2 (69%) optical characteristics are E-like. 
This preponderance of E-like hosts among AGN mergers is also 
seen when using the ISO classification (Fig. 13£>; 8 out of 10 or 
80% of ISO AGNs are E-like). A similar trend is seen among 
warm objects (10 out of the 13 objects with /2s/ ho > 0.2 or 
77% are E-like; Fig. 13c). Starburst and LINER systems, on 
the other hand, show no preference between disk-like and E- 
like hosts. 

The nature of the host does not seem to depend strongly 
on the infrared luminosity, the 60 /xm/100 /im ratio, or the 
strengths of the stellar H/3 and Mg lb features. But, as expected 
in the merger scenario for elliptical galaxy formation, E- and 
E/D profiles are more common in the later phases of interac- 
tion [class V (27%) and especially class IVb (49%); Fig. 13a"], 
while ambiguous profiles are more often seen in early mergers 
(IVa, 59%), a sign perhaps that many of these systems are still 
showing the damages caused by the recent merger event. 

4.2. Absolute Magnitudes 

The absolute R and K' magnitudes of the host galaxies in the 
1-Jy sample are listed in the last two columns of Table 2. Me- 
dian and mean values are listed in Table 3. The distributions 
for double- and single-nucleus systems are shown separately in 
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Figure 14. The median and mean /?-band luminosities of the 
individual galaxies involved in binaries are 0.86 and 0.85 ± Jj 43 
L* (1 a), respectively. At K', the median and mean for these ob- 
jects are 0.91 and 0.90 ± J L*. For single-nucleus systems, 
the median and mean /?-band luminosities are 1.69 and 1.69 
± L*, and the median and mean /T'-band luminosities are 
1.92 and 2.36 ± \ Previous studies have obtained similar 
host galaxy magnitudes: <~ 2.5 L B (Armus, Heckman, & Miley 
1990), - 3 L* (ranging from ~ 1 - 24 L*) for the ULIGs in the 
2-Jy sample (Murphy et al. 1996), ~ 2.4 L* H for warm ULIGs 
from the IRAS BGS sample (Surace & Sanders 1999), and - 1 - 
2.5 L* H for cool ULIGs from the BGS and 1-Jy samples (Surace, 
Sanders, & Evans 2000). Colina et al. (2001) derive noticeably 
fainter, sub-L*, host galaxy magnitudes (< Mh > = -23.5 ± 
0.7 mag. scaled down to our cosmology) for 27 mostly cool 
ULIGs. The //-band HST/WFPC2 data used for this analysis 
are shown in Bushouse et al. (2002). Direct comparisons with 
our data for the four objects in common with the 1-Jy sample 
indicate that much of the diffuse low surface brightness features 
in our data are not visible in the snapshot images of Bushouse 
et al. (2002). This underestimate of the light contribution from 
the host galaxies explains the fainter host galaxy magnitudes 
derived by Colina et al. (2001). 

The large number of galaxies in the 1-Jy sample allows us to 
look for trends between host galaxy luminosities and infrared 
luminosities, optical spectral types, and IRAS colors. Focussing 
on single-nucleus systems, there appears to be a tendency to 
find luminous hosts in the most infrared luminous systems of 
the 1-Jy sample (see Fig. 15a and Table 3). This is due to the 
fact noted in §2.1 that systems with L 1R > 10 12,3 L are more 
distant on average (< z > = 0.183 ± 0.053) than ULIGs with 
L m < 10 12 ' 3 L Q (< z > = 0.132 ± 0.029). 

Figure 15 and Table 3 indicate that the most luminous hosts 
also tend to be found in systems with warm IRAS 25 yum/60 fim 
colors (> 0.2) and AGN characteristics. The hosts of the warm 
objects are 0.6 - 0.7 magnitudes brighter on average than the 
hosts of the cool objects, although both distributions overlap 
considerably. In this case, the difference in mean host lumi- 
nosity cannot be attributed to a difference in redshift since both 
classes of objects have virtually the same median redshift (= 
0.148). The same is true when comparing the Seyfert and non- 
Seyfert populations (< z > = 0.143 ± 0.056 and 0.150 ± 0.041, 
respectively), yet the host galaxies of Seyferts are brighter by 
about 0.8 - 1 .0 mag. on average at both R and K' than the hosts 
of LINERs and H II region-like galaxies. K-S tests between 
the distributions of Seyferts vs non-Seyferts and warm vs cool 
ULIGs confirm these differences (P[null] < 10" 5 and 10" 2 , re- 
spectively). These results can be attributed to real differences 
in the host luminosities of AGNs vs non-AGNs, or may be due 
to a prominent starbursting population in warm Seyferts con- 
tributing significantly at R and K' . 

One should be cautious when interpreting these results, how- 
ever. An object-by-object comparison with the //-band data 
of Surace & Sanders (1999) suggests that the host luminosi- 
ties measured in a few of the warm Seyferts appear to be 
slightly higher than those measured from AO imaging (assum- 
ing H-K' < 0.5 for the host; see below). This may indicate 
that we are slightly underestimating the contribution from the 
central point source in some of these objects (the possibility 
that the AO imaging of Surace & Sanders is underestimating 
the host luminosities appears less likely because the exposure 
times for these observations were fairly long and therefore the 



surface brightness limits of the data were fairly deep). The PSF 
subtraction in our data is more difficult in the warm, highly 
nucleated systems of the 1-Jy sample than in the cooler, more 
diffuse systems; this may introduce a slight bias which may 
account for the weak trends between host luminosities, IRAS 
colors and spectral types. These issues are discussed in more 
detail in §§4.5, 5.2, and 5.3. 

Next we compare the host magnitudes of ULIGs with the 
results from studies of quasar host galaxies. At R, the host 
galaxies in ULIGs appear slightly less luminous on average 
than the HST measurements of Dunlop et al. (2002) on radio- 
quiet quasars (-22.35 ± 0.15 mag., adjusted to H = 75 km s" 1 
Mpc" 1 ) and radio-loud quasars (-22.85 ± 0.10 mag.). Note, 
however, that ULIG hosts span a broad range in luminosity 
so there is considerable overlap between ULIGs and the other 
classes of objects. Given the correlations between host lumi- 
nosities, infrared luminosities, IRAS 25 yum/60 /im colors, and 
optical spectral types discussed in the previous paragraphs, we 
find that if we focus our analysis on the subset of objects with 
log [Li R /L ] > 12.3, /25 / /60 > 0.2, and with Seyfert character- 
istics at optical wavelengths, the apparent shift in host galaxy 
R magnitudes between these ULIGs and the quasar dataset of 
Dunlop et al. disappears. 

To directly compare our data with the ground-based results 
of McLeod & Rieke (1994a, 1994b) and HST/NICMOS data of 
McLeod & McLeod (2001) on quasars, we need to transform 
the //-band magnitudes listed in these papers into K' magni- 
tudes. For an old stellar population at z = 0, H -K' = Q.\6 
(Surace 1998), but H-K 1 w 0.3-0.4 if < z >= 0.1-0.2 (e.g., 
Lilly & Longair 1984). Given these corrections, the low- 
luminosity sample of quasars of McLeod & Rieke (1994a) have 
an average K' host magnitude Mk> ~ -24.3 ±0.6, while the 
high-luminosity quasars have Mk' ~ -25.1 ±0.7. The host 
galaxies in the 1-Jy ULIGs are therefore comparable in lumi- 
nosity to those of the /z/g/z-luminosity quasars. So even though 
we found in §2.1 that the 1-Jy ULIGs have lower integrated 
K' luminosities than the high-luminosity quasars of McLeod & 
Rieke (1994b), we find that they have host galaxies with sim- 
ilar K' luminosities. This good agreement is also found in the 
more recent analysis of the high-luminosity sample by McLeod 
& McLeod (2001). 

The claim by McLeod & Rieke (1994a, 1994b; ground- 
based //-band data) and McLeod, Rieke, & Storrie-Lombardi 
( 1 999; HST/NICMOS data) of a trend between nuclear and host 
galaxy luminosities in quasars has recently come under ques- 
tion. The /?-band results of Dunlop et al. (2002) do not seem 
to confirm this trend. A recent Zf-band ground-based study of 
14 very high luminosity {My < -25) radio-quiet quasars with 
0.26 > z > 0.46 by the same group (Percival et al. 2001) 
find < M K >« -24.3 ± 0.04, a value closer to the host galaxy 
magnitudes of the tow-luminosity quasar sample of McLeod 
& McLeod (1994a) and significantly less luminous on average 
than the host galaxies of the 1-Jy ULIGs. However, as stated 
for instance by Dunlop et al. (2002) when discussing the study 
of Taylor et al. (1996), the results from ground-based studies 
of quasars beyond z = 0.15 like that of Percival et al. should be 
treated with caution given the technical difficulties involved in 
removing the nuclear PSF to characterize the underlying host 
galaxy of these objects (see also discussion in §4.5). 

Adaptive optics imaging of quasars alleviates the problems 
associated with PSF subtraction in ground-based data. The re- 
cent study by Surace et al. (2001) of 18 PG quasars with strong 
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far-infrared excesses finds //-band host galaxy luminosities that 
range from 0.5 to 7.5 L* , with a mean of 2.3 L* , consistent with 
the average K 1 host galaxy luminosity of 1-Jy ULIGs measured 
in the present study and similar to that found by McLeod & 
McLeod (2001) among quasars. 

4.3. Luminosity Ratios in Binary Systems 

The mass ratio of a binary system of galaxies is an impor- 
tant parameter in predicting the amount of "damage" incurred 
in the interaction. Minor mergers comprised of objects with 
mass ratios of order ~ 10 are predicted to have only benign 
effects on the main component of the system (e.g., Walker, Mi- 
hos, & Hernquist 1996; Bendo & Barnes 2000; Naab & Burkert 
2001). Major mergers with mass ratio of ~ 3 or less appear to 
be needed to create a ULIG (e.g., Mihos & Hernquist 1996). 
The kinematic properties of ULIGs also appear consistent with 
this scenario (e.g., Genzel et al. 2001). Given our lack of de- 
tailed kinematic data, we have no reliable way to determine the 
total mass of each galaxy in our binary systems. We make the 
assumptions that the total masses scale linearly with the stellar 
masses and that the stellar masses scale with the R or K' lumi- 
nosities of the host galaxies. Figure 16 shows the distributions 
of AM(R) and AM(K') (equivalent to luminosity ratios) of the 
Class III systems in the 1-Jy sample for which we were able 
to reliably separate each galaxy (i.e. most of them are Class 
Ilia). Most systems have a magnitude difference of less than 1 .5 
or a luminosity ratio of less than 4, in general agreement with 
the predictions from numerical simulations of galaxy mergers 
if mass indeed scales linearly with luminosity. The magnitude 
difference is slightly larger at K' than at R. The two objects 
with the largest K' magnitude differences are F00456-2904 and 
F12127-1412. 

4.4. Colors 

In Figure 17, we find no significant correlation between 
the nuclear color and the underlying galaxy color (R - K') host 
(P[null] > 0.01 for H II, LINER, and Seyfert galaxies). In con- 
trast to the nuclear colors, the median color of the underlying 
galaxies of ULIGs is only slightly redder than that of normal 
galaxies or the host galaxies of quasars [(R-K')^^ = 3.0 for 
either H II galaxies, LINERs, or Seyfert 2 galaxies compared 
to 2.6 and 2.5 for normal galaxies and quasars, respectively; 
de Vaucouleurs & Longo 1988; Dunlop et al. 2002]. This re- 
sult confirms that the obscuring gas and dust is preferentially 
concentrated within the inner 4 kpc of ULIGs (e.g., Scoville 
et al. 2000). The slightly redder colors of Seyfert 1 hosts 
[with median (/?-/f')host of 3.4] may be due to contamination 
at K' by hot dust emission from the AGN or residuals from 
the bright central source subtraction procedure. The colors of 
ULIG hosts are similar to those found by Surace et al. (2001) 
among infrared-excess PG quasars. 

4.5. Half-Light Radii and Surface Brightnesses 

The seeing-deconvolved half-light radii, r e , of the elliptical- 
like host galaxies in our sample are listed in Table 2 and plot- 
ted in Figure 18. The median and mean of r e at R are 4.60 
kpc and 4.80 ± 1.37 kpc (1 a), while at K' we find 3.55 kpc 
and 3.48 ± 1.39 kpc (1 a). These half-light radii are virtu- 
ally indistinguishable from those measured at H by McLeod & 
McLeod (2001) in the eight E-like "high-luminosity" quasars 
of their sample (< r e > = 3.39 ± 1.90 kpc; Fig. lSb). On the 
other hand, the values measured by Dunlop et al. (2002) among 



high-luminosity radio-quiet quasars (7.63 ±1.11 kpc adjusted 
to Hq = 75 km s -1 Mpc -1 ) and radio-loud quasars (7.82 ± 0.71 
kpc) appear to be significantly larger than those measured in 
the 1-Jy ULIGs, although there is some overlap between the 
two distributions. We do not try to compare our results with 
the half-light radii measured by Taylor et al. (1996) at K be- 
cause these data have proven to be unreliable: An object-by- 
object comparison between the data of Taylor et al. (1996) and 
Dunlop et al. (2002) reveal discrepancies of up to a factor of 
nearly <~ 25 (e.g., 1549+203) in r e , which cannot be attributed 
to color gradients (R-K 1 ) but rather are probably due to incor- 
rect PSF subtraction in the ground-based data of Taylor et al. 
Similar technical difficulties are likely to affect to some extent 
the ground-based data of Percival et al. (2001), given the rela- 
tively high redshifts of their sources. 

Combining measurements of the half-light radius and mean 
/?-band surface brightness within half-light radius, [i e , Dunlop 
et al. (2002) have argued that the host galaxies of luminous 
quasars follow the same fi e - r e relation as that of normal (inac- 
tive) ellipticals. The surface brightnesses of ULIGs listed in Ta- 
ble 2 have been corrected for cosmological dimming [oc (1 +z) 4 ] 
and are compared with those of elliptical galaxies and quasars 
in Figure 19a. The elliptical-like (E and E/D) host galaxies of 
the 1-Jy ULIGs have slightly brighter /?-band surface bright- 
nesses than the hosts of the quasars studied by Dunlop et al. 
Nevertheless, the elliptical-like host galaxies of ULIGs appear 
to follow the n e - r e relation of normal ellipticals, although with 
considerable scatter. In this figure, we used the fit of Hamabe 
& Kormendy (1989) and shifted it in magnitude by V-R = 0.5 
mag (keeping the slope constant). Seyfert 1 galaxies show a 
tendency to fall slightly above the data for non-Seyfert Is. The 
source of this tendency is discussed below. 

Interpretation of our data at K' is more difficult because r e 
and yi e have been measured reliably in only 18 ULIGs of our 
sample (this analysis was carried out only when the K' images 
were deep enough to reliably measure the surface brightness 
profiles out to a radius comparable to the optical radius). The 
results are presented in Figure 19b. A weak anti-correlation is 
found between apparent surface brightnesses and galaxy sizes. 
This anti-correlation is stronger if we exclude Seyfert 1 galaxy 
hosts from the fit (solid line in Fig. 19b). The data for the 
Seyfert Is should be considered with caution since they are sig- 
nificantly more sensitive to our PSF subtraction procedure than 
for the other objects. The bright [i e of the three Seyfert 1 host 
galaxies near the top of Figure 19b may be due to positive resid- 
uals from the AGN. 

We first compare our K 1 results with the V-band [i e - r e rela- 
tion of Hamabe & Kormendy (1987) for normal ellipticals and 
find a generally good agreement with the non-Seyfert 1 data 
if we assume a constant color V —K' = 3.7 mag. This color 
term is slightly redder than that of normal elliptical galaxies 
(e.g., Nolan et al. 2001; Fukugita et al. 1995), possibly be- 
cause the hosts of ULIGs are dustier than elliptical galaxies and 
star formation (red supergiants) contributes at K' . [We did not 
use the B-band data of Bender, Burstein, & Faber (1992) and 
Faber et al. (1997) for this comparison to avoid the even larger 
color correction.] We also compare our results with the ZT-band 
data of Pahre (1999). The slope of the yi e -r e relation derived 
from these data (dotted line in Fig. 19V) is shallower than that 
of Hamabe & Kormendy (1987). Unfortunately, a comparison 
with quasar hosts is not possible because of the lack of reliable 
quasar host data in the near-infrared (McLeod & McLeod 2001 
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list r e but not [i e in their paper, and the data of Taylor et al. 
are not considered reliable for the reasons discussed in the first 
paragraph of this section). 

There is no obvious correlation of the /?-band r e (or fj, e ) 
with IRAS 25 /im/60 fim and 60 /im/100 fim colors, infrared 
luminosities, optical and mid-infrared spectral types, and the 
strengths of the stellar H/3 and Mg lb absorption features. There 
are too few measurements of r e at K' to search for statistically 
significant trends with this parameter. Recall, however, the ex- 
istence of correlations between the host galaxy R and K' lumi- 
nosities (oc fi e r^) and infrared luminosities, IRAS 25 yum/60 fim 
colors, and AGN characteristics noted in §4.2 (Fig. 15 and Ta- 
ble 3). 

4.6. Bars, Axial Ratios, and Boxiness 

The ~ 1" spatial resolution of our images prevents us from 
studying in detail the cores of ULIGs on scales smaller than a 
few kpc. However, a few statements can be made on galactic 
structures seen beyond the seeing disk. Only three objects in 
our sample, F14394+5332 W (triplet), F23233+2817 (isolated 
spiral), and F23327+2913 S (widely separated pair), show ev- 
idence for bar-like structures which extend 8 kpc (P.A. = 90°), 
12 kpc (20°), and 11 kpc (120°), respectively. None of the 
single-nucleus merger remnants shows a strong large-scale bar. 

Our ellipse-fitting routine of the surface brightness distribu- 
tion provides quantitative information on the axial ratio (b/a) 
of the host galaxies on scales larger than twice the seeing disk. 
Figure 20 shows the radial dependence of the axial ratio for 
the 58 single-nucleus mergers of our sample with "E" or "E/D" 
profiles (F23233+2817 is not included in this subset of objects 
because it is an isolated spiral galaxy). In the majority of cases, 
the axial ratio is constant with radius to within ±0.15. We plot 
in Figure 21 the distribution of b/a measured at 3 x r e . The 
median and mean of this quantity are 0.77 and 0.75 ± 0.15 (1 
a), respectively. This is similar to the average values measured 
in normal elliptical galaxies (e.g., Ryden 1992) and in the hosts 
of quasars (Dunlop et al. 2002). 

Surface photometry of elliptical galaxies often shows 
isophotes that deviate slightly from a perfect ellipse. The depar- 
ture of the isophotes from a perfect ellipse is measured by the 
A\ coefficient of the cos 4<fi term in equation (1). Typical val- 
ues for ellipticals normalized to the major-axis length, a, range 
from -0.02 to +0.03 (e.g., Bender et al. 1989), where the neg- 
ative values indicate boxy isophotes and positive values indi- 
cate disky isophotes. Ellipticals with boxy isophotes are gener- 
ally luminous (massive) radio-loud ellipticals with X-ray halos, 
very slow stellar rotation and a distinct cuspy cores, whereas 
ellipticals with disky isophotes are low-to-moderate luminos- 
ity (mass) lenticular and elliptical galaxies which are generally 
radio-quiet, show little or no X-ray emission, and exhibit sig- 
nificant rotational support and a power-law surface brightness 
profile on small scales (e.g., Kormendy & Bender 1996; Faber 
et al. 1997; Pahre, de Carvalho & Djorgovski 1998 and refer- 
ences therein). 

The radial profiles of the normalized /?-band boxiness pa- 
rameter in the single-nucleus ULIGs of our sample are shown 
in Figure 22. Several of these galaxies have large absolute val- 
ues of A\ja (generally from -0.05 to +0.05) compared to nor- 
mal ellipticals. Significant A\ja gradients are seen in many ob- 
jects. Thus it is hard to attribute a single boxiness value to each 
galaxy. The distributions of An,/ a measured at 2 x r e and 3xr f 
in elliptical-like (E and E/D) ULIGs are presented in Figure 23. 



There is no clear preference between positive and negative box- 
iness values. The median and mean of An/a at 2 x r e are -0.002 
and -0.001 ± 0.036, and +0.005 and +0.014 ± 0.059 at 3 x r e . 
The broad range of boxiness values in ULIG hosts probably re- 
flects the residual effects of the violent merger that took place 
in these galaxies. 

5. DISCUSSION 

In this section, the results from our imaging survey are 
combined with those from other published studies on ULIGs, 
quasars, and normal elliptical galaxies to attempt to answer the 
following questions: (1) What is the nature of ultraluminous 
infrared galaxies?; (2) Are ultraluminous infrared mergers el- 
liptical galaxies in formation? (3) Are ultraluminous infrared 
galaxies quasars in formation? 

5.1. What is the Nature of Ultraluminous Infrared Galaxies ? 

What types of galaxies make up ultraluminous infrared sys- 
tems? What conditions are needed to trigger their quasar-like 
infrared luminosities? Answers to these two questions would 
go a long way in clarifying the nature of ULIGs. There is now 
a consensus among researchers in this field concerning the fre- 
quency of interactions in ULIGs and the number of galaxies in- 
volved in these interactions. The results from the present survey 
(§3.1) show convincingly that the great majority of ULIGs (~ 
95%) are involved in the interaction of two, and no more than 
two, galaxies. The recent HST results of Bushouse et al. (2002) 
also support this conclusion. The higher frequency of multiple 
mergers claimed earlier by Borne et al. (2001) and Cui et al. 

(2001) has now been shown to be due primarily to the presence 
of several star-forming super star clusters which mimic low- 
luminosity nuclei at optical wavelengths. Near-infrared obser- 
vations reveal the true nature of these clusters, ruling out gen- 
uine galactic nuclei. Follow-up spectroscopy of apparent neigh- 
bors to ULIG systems also reveals in some cases that they are 
not physically associated with the ULIGs (see Appendix A). 

Our detailed morphological analysis of the 1-Jy sample pro- 
vides new statistical evidence that most (56%) ULIGs are in- 
volved in the late phase of a merger where the two nuclei have 
merged into a single object. This frequency does not seem to 
depend strongly on spatial resolution since the HST data of 
Scoville et al. (2000), Colina et al. (2001), and Bushouse et al. 

(2002) reveal very few double-nucleus systems that were not al- 
ready known from arcsecond-resolution ground-based data like 
ours. As described in §3.2, the frequency of advanced merg- 
ers is a function of infrared luminosity, increasing dramatically 
above ~ 10 12,5 Lq and decreasing among luminous infrared 
galaxies with log [L IR /L Q ] = 1 1 .50 - 1 1 .99. So not only can we 
say that the ultraluminous infrared phase is generally triggered 
when the individual galaxies in the merger have lost their indi- 
vidual identity and the nuclei have fused into a single entity, but 
also that the extreme environment during the final phases of a 
merger appears to be a necessary condition to produce the most 
extreme examples of ULIGs. 

Another ~ 39% of the ULIGs in the 1-Jy sample harbor two 
distinct nuclei but always show tidal tails at R, evidence that 
the ultraluminous infrared phase here was triggered after the 
first close passage of the two interacting galaxies. These re- 
sults are at least qualitatively consistent with the simulations 
of Mihos & Hernquist (1996) where two merging galaxies go 
through two distinct episodes of enhanced star formation ac- 
tivity (and therefore large infrared luminosity; note that Mihos 
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& Hernquist do not try to model AGN activity) with the first 
peak taking place shortly (~ 1 - 2 x 10 8 yr) after the first 
close encounter, when the two galaxies are still widely sepa- 
rated, and the second peak taking place during the final merger 
of the galaxies, ~ 10 9 yr after the first encounter. Our results 
indicate that in quite a few cases the star formation activity trig- 
gered during the first passage is sufficient to reach the ULIG 
threshold (although we have no high-resolution maps in the in- 
frared to determine whether the FIR-emitting region in these 
widely spaced ULIG systems is a kpc-scale ring, as predicted 
by the simulations of Mihos & Hernquist). The detailed kine- 
matic study of four ULIGs (including F10190+1322 from the 
1-Jy sample) by Murphy et al. (2001a) confirms that the ULIG 
phase may sometimes take place well before the final merger. 

In the simulations of Mihos & Hernquist, the relative am- 
plitude of the peaks of enhanced star formation activity is 
only weakly dependent on the initial orientation of the merg- 
ing galaxies but depends strongly on the morphological proper- 
ties of these galaxies, particularly their bulge-to-disk mass ratio. 
Unfortunately, reliable profile decomposition to determine this 
ratio in our sample of binary galaxies requires sub-kpc resolu- 
tion and therefore cannot be done with the present data. Nev- 
ertheless, our data analysis provides some constraints on the 
type of galaxies involved in the interactions. Our crude pro- 
file classification in §4. 1 suggests that an important fraction of 
these binary galaxies have a prominent spheroidal component 
(possibly a bulge). As discussed in §4.2 (and in more detail 
in §5.2 and §5.3 below), the host galaxies of single-nucleus 
ULIGs span a broad range in luminosities with average of or- 
der ~ 2 L* . Binary systems involve two <~ 0.9 L* (on average) 
galaxies with luminosity ratio generally less than 4. This is ev- 
idence that mergers between small galaxies or minor mergers 
involving the accretion of small satellite galaxies by larger ones 
do not provide the conditions necessary to trigger an ultralumi- 
nous infrared phase. Large torques appear needed to produce 
enough "damage" to the galaxies involved in the interaction so 
that gas transfer to the inner portions of the galaxies is sufficient 
to trigger a strong episode of star formation and/or AGN activ- 
ity and reach the ULIG phase. These results are qualitatively 
consistent with numerical simulations of unequal mass mergers 
(e.g., Mihos & Hernquist 1996; Walker, Mihos, & Hernquist 
1996; Bendo & Barnes 2000; Naab & Burkert 2001) and the 
recent kinematic data of Genzel et al. (2001). 

There are interesting exceptions to the merger scenario for 
ULIGs. Among these is F23233+2817, the only spiral galaxy 
in the 1-Jy sample which does not seem to be involved in any 
galactic interaction. The infrared luminosity of this object lies 
on the boundary between luminous and ultraluminous infrared 
galaxies (log [Li R /L Q ] = 12.00), so it may be that the Seyfert 2 
nucleus in this object (or unsuspected circumnuclear starburst) 
is sufficiently powerful to heat enough dust to reach the ultra- 
luminous threshold without the need for refueling from a major 
merger. Eight ULIGs in the 1-Jy sample appear to involve inter- 
acting galaxies that are widely (> 20 kpc) separated (see Table 

1 and Fig. 5; one of these objects, Fl 1223-1244, hosts a Seyfert 

2 nucleus). All of these objects show tidal features indicative 
of a pre-merger phase (i.e. after first contact; this seems to rule 
out the possibility that these objects are observed before first 
contact, as suggested by Dinh-V-Trung et al. 2001). It would 
therefore appear that an ultraluminous infrared phase can also 
be triggered when two galaxies in a pair are near apgalacticon 
(farthest approach after first passage). This is difficult to recon- 



cile with the predictions of the merger scenario. A way out of 
this uncomfortable situation is to suppose that the ULIG phase 
was triggered by an earlier merger event which is not evident at 
the present epoch. One should also be cautious when interpret- 
ing the data on these eight ULIGs since the physical connection 
between the galaxies in most of these pairs has not yet been 
confirmed through detailed spectroscopy (see Appendix A). 

5.2. Are Ultraluminous Infrared Mergers Elliptical Galaxies 
in Formation? 

Among the 66 single-nucleus mergers in the 1-Jy sample, 48 
(73%) of them are fit adequately by a de Vaucouleurs profile at 
R and K' (§4.1). This fraction is uncertain, however, since the 
relatively modest dynamical range of our data (especially at K') 
prevents us from deciding if a de Vaucouleurs profile is a better 
fit than an exponential disk in 25 (38%) of these objects (these 
are the so-called E/D galaxies in Table 2 and Fig. 12). Nev- 
ertheless, we can state with confidence that elliptical-like hosts 
are observed in the majority (> 50%) of late mergers of the 
1-Jy sample. This result supports the conclusions from previ- 
ous studies that were based on smaller samples of ULIGs (e.g., 
Wright et al. 1990; Stanford & Bushouse 1991; Doyon et al. 
1994; Scoville et al. 2000; Cui et al. 2001), and is at least qual- 
itatively consistent with the scenario that ULIGs evolve into 
ellipticals through merger-induced, dissipative collapse (Kor- 
mendy & Sanders 1992). In the rest of this section, we further 
explore the predictions of this scenario by comparing the prop- 
erties of the elliptical-like hosts found in the 1-Jy sample with 
those of nearby spiral bulges, lenticular (SO) galaxies, and el- 
lipticals, while taking into account the various sources of errors 
in our analysis of the ULIGs. 

In §4.2, the elliptical-like (E and E/D) hosts in the 1-Jy sam- 
ple were found to have a median and average luminosities at 
R of 1.77 and 1.66 ± L*, and 2.00 and 1.98 ± J;| L* at 
K' . Residuals from AGNs and compact circumnuclear star- 
bursts may affect (overestimate) slightly the luminosities of 
some of these hosts - this effect appears to be more important 
at K' than at R and more frequent among Seyfert Is than among 
Seyfert 2s (§4.2 and §4.5). The results at R are thus consid- 
ered more reliable, even though the effects associated with dust 
extinction may be more significant here than at near-infrared 
wavelengths. Assuming V-R ~ 0.6 mag for typical ellipticals 
(Fukujita et al. 1995), the average V-band absolute magnitude 
of the elliptical-like hosts in the 1-Jy sample is ~ -21.1 ± 0.7 
mag. The elliptical-like hosts of ULIGs therefore span the in- 
termediate luminosity range between the low-luminosity (< 1 .0 
L*) "power-law galaxies" and high-luminosity (> 3.5 L*) "core 
galaxies" of Faber et al. (1997). 

The distribution of axial ratios among elliptical-like hosts in 
the 1-Jy sample resembles that of elliptical galaxies (§4.6). The 
surface brightness profiles of ULIG hosts also follow the /i c - 
r e projection of the fundamental plane of nearby ellipticals, but 
considerable scatter is present (§4.5; Fig. 19). The average r e 
of the 1-Jy ULIGs (4.80 ± 1.37 kpc at R and 3.48 ± 1.39 kpc 
at K') are fairly typical of intermediate-luminosity ellipticals 
(e.g., Kormendy & Djorgovski 1989; Bender et al. 1992; Faber 
et al. 1997). We mentioned in §4.5 that residuals from the PSF 
subtraction may artificially increase the values of the average 
surface brightness within half-light radius (the half-light radius 
may also be slightly Mntferestimated). This effect is likely to 
be at the origin of some of the scatter in Figure 19 and may 
also explain why some Seyfert 1 hosts in our sample lie signif- 
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icantly above the [i e -r e relation of ellipticals. Dust extinction 
and emission and circumnuclear star formation are other likely 
sources of scatter in this diagram. 

Given the recent violent merger events in ULIGs, it is per- 
haps surprising to find that any of these galaxies follow the 
/Lt c - r e of ellipticals. This result brings support to numerical 
simulations which show that violent relaxation in dissipative 
gas-rich merger systems is very efficient. A general rule of 
thumb pointed out by Mihos (1999) is that "violent relaxation 
and/or dynamical mixing occurs on a few rotation periods at the 
radius in question?' This time scale is only ~ 10 8 yrs in the cen- 
tral (< kpc) regions of the mergers, but exceeds 10 9 yrs in the 
outer (> a few r e ) regions. Indeed, while the central regions of 
many ULIG hosts appear to mimic those of elliptical galaxies, 
the effects of the mergers still appear to be visible beyond a few 
r e as illustrated by the broad range of boxiness values measured 
in ULIG elliptical-like hosts of the 1-Jy sample (§4.6; Fig. 23). 
The large fraction (26%) of ULIG hosts which are not fit ad- 
equately by either an exponential or de Vaucouleurs profile is 
another indication that the light distribution in several objects 
may still be showing the effects of the recent mergers (patchy 
dust distribution could also explain these results). 

It will be important in the future to obtain deep high- 
resolution (< 100 pc) images of the 1-Jy ULIGs in both the 
visible and near-infrared domains to improve on our AGN/PSF 
subtraction and determine if the ULIG host galaxies are power- 
law or core galaxies with properties that are consistent with the 
photometric portion {ji e -r e ) of the core and global fundamental 
planes of early-type galaxies (e.g., Faber et al. 1997). Multi- 
color maps combined with detailed radiative transfer models 
will be particularly useful in these objects to help disentangle 
the effects of dust extinction and emission and circumnuclear 
star formation on this scale (see, for instance, Regan, Vogel, & 
Teuben 1995). ULIGs are known to have sub-kpc nuclear con- 
centrations of molecular gas with mean surface densities sim- 
ilar to the stellar densities measured in the cores of ellipticals 
(a few x 10 4 M Q pc" 2 ; e.g., Bryant & Scoville 1999; Downes 
& Solomon 1998). If post-VLlG hosts are indeed genuine el- 
lipticals, one might therefore expect the hosts of current star- 
forming ULIGs to have lower stellar densities than the core 
densities of classical ellipticals. This effect could be signifi- 
cant since the star formation rates of ULIGs are <~ few x 100 
M yr" 1 and the ULIG phase may last 10 7 - 10 8 years. The 
effect at large radii should be less important because most of 
the star formation (and AGN) activity in ULIGs is concentrated 
well within the central kpc (e.g., Soifer et al. 2000, 2001). 

Given that merger systems are expected to reach their equi- 
librium values of light distribution and kinematics on a short 
time scale (e.g., Mihos 1999; Bendo & Barnes 2000), the kine- 
matic properties of ULIGs should also follow the fundamental 
plane of ellipticals. Considerable efforts have been invested in 
recent years to constrain the stellar kinematics in the cores of 
luminous infrared mergers (e.g., Doyon et al. 1994; Shier & 
Fisher 1998; James et al. 1999; Genzel et al. 2001). The ad- 
vent of powerful new near-infrared spectrographs on 8-meter 
class telescopes now allows the radial profiles of the velocity 
dispersion and rotation velocity to be measured in the brightest 
systems and compared with those of elliptical galaxies. Most 
relevant to our study of the 1-Jy ULIGs are the results from 
Genzel et al. (2001) since they refer to ULIGs while those of 
Shier & Fisher (1998) and James et al. (1999) refer to systems 
of lower luminosities. Genzel et al. (2001) find that random 



motions dominate the stellar dynamics in the cores of their (12) 
systems, but that significant rotation is also common. The prop- 
erties of these galaxies closely resemble those of intermediate- 
mass (<~ L*), disky ellipticals. These results are consistent with 
ours in the fact that they rule out giant (> 3.5 L*) ellipticals 
with large cores and little rotation as the hosts of ULIGs. A 
more quantitative comparison is difficult due the broad range 
of photometric properties found in the 1-Jy sample and the rel- 
atively small size of the kinematic sample of Genzel et al. Re- 
cently commissioned three-dimensional spectrographs on large 
telescopes should provide the sensitivity and two-dimensional 
spatial coverage needed to derive detailed kinematic maps for 
several more ULIG hosts, and therefore allow more robust com- 
parisons with our results. 

5.3. Are Ultraluminous Infrared Galaxies Quasars in 
Formation? 

The results presented in §2 and §3 indicate clear trends be- 
tween merger phase, IRAS colors, and the importance of AGN 
activity in ULIGs (in addition to the trend between infrared 
luminosity and merger phase, which has already been dis- 
cussed in §5.1). Objects with warm quasar-like infrared col- 
ors (as measured by the IRAS 25 /im/60 /im flux ratio) present 
strong signs of AGN-dominated activity at both optical and 
mid-infrared (ISO) wavelengths and are generally found in ad- 
vanced mergers (based not only on the apparent nuclear separa- 
tion but also on the overall morphology of the system including 
the shape and compactness of the nuclear core and the strengths 
and lengths of the tidal tails). These trends are consistent with 
an evolutionary sequence in which cool starburst-dominated 
ULIGs transform into warm AGN-dominated ULIGs as the 
merger of the two <~ 0.9 L* galaxies takes place. 

It is clear that exceptions to these general trends with merger 
phase abound in the 1-Jy sample. For instance, 46% of the 41 
advanced mergers (classes IVb and V in §3.4) in our sample 
show no obvious signs of Seyfert activity (see Fig. 10). There 
are also seven Seyfert 2 nuclei (but no Seyfert Is) among the 
45 pre-mergers (classes Ilia and Illb). The cases of an isolated 
AGN with a spiral galaxy host and a widely separated pair with 
an AGN were mentioned in §5.1. Although merger phase ap- 
pears to be an important factor in determining the relative im- 
portance of star formation and AGN activity in ULIGs, it is not 
the only one (see the discussion in Farrah et al. 2001). Star 
formation activity depends on the molecular gas supply and the 
efficiency to transform this gas into stars (which is a compli- 
cated function of density, metallicity, magnetic filed strength, 
local gas kinematics, etc.). On the other hand, the level of AGN 
activity in ULIGs depends linearly on the local mass accretion 
rate onto the supermassive black holes that are presumed to ex- 
ist in these objects, and the efficiency of the system to transform 
the gravitational energy into radiation that can escape the sys- 
tem and be detected by our instruments. One only needs an 
accretion rate of ~ 1 M yr" 1 to produce the infrared (bolo- 
metric) luminosity of a ULIG (assuming a radiative efficiency 
of 10% in rest-mass units). Given the complexity of the star 
formation and AGN phenomena it is not surprising to see large 
scatter in the correlations between merger phase and AGN ac- 
tivity. Note, however, that a late merger phase appears to be 
a necessary condition to trigger powerful Seyfert 1 (= quasar) 
activity. Our results emphasize the importance of studying a 
large and homogeneous set of ULIGs like the 1-Jy sample to 
draw statistically meaningful conclusions on these objects. The 
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problems of small sample sizes and/or inhomogeneous selec- 
tion criteria have plagued many studies of luminous infrared 
galaxies in the past. 

The natural end-result of the merger sequence described in 
the first paragraph is an object with quasar-like colors and 
luminosities that shows only slight residual effects from the 
merger. This object has several properties in common with 
nearby quasars. A summary of the properties of quasar host 
galaxies is beyond the scope of this paper. Suffice it to say that 
evidence for galaxy interactions and mergers is seen in some 
but not all nearby quasars (see, e.g., Hutchings & Neff 1992; 
Hutchings et al. 1994; Disney et al. 1995; Bahcall et al. 1997; 
Stockton 1999; Dunlop et al. 2002 and references therein). The 
present discussion focusses on the so-called infrared-excess 
quasars since they are particularly good candidates for transi- 
tion objects which have just gone through a ULIG phase and 
are now in the process of settling in to become optical quasars. 
These objects show morphological evidence for recent merg- 
ers (e.g., Surace et al. 2001 and references therein), contain 
large quantities of molecular gas (e.g., Evans et al. 2001 and 
references therein), and sometimes present spectroscopic evi- 
dence for substantial gas outflows (e.g., BAL phenomenon, Li- 
pari, Colina, & Macchetto 1994). The host galaxies of these 
infrared-excess quasars have luminosities and colors which are 
similar to those of the 1-Jy ULIGs (§4.2 and §4.3). In a series of 
papers, Canalizo & Stockton (2000a, 2000b, 2001) have shown 
that these "transition" quasars show strong recent star-forming 
activity, which in most case is directly related to the tidal in- 
teraction. The ages they derive for the starburst populations 
range from currently active star formation in some objects to 
poststarburst ages <~ 300 Myr in others. They argue for a direct 
connection between interactions, starbursts, and QSO activity. 

The issue of the hosts in ULIGs and quasars is important 
and deserves further attention. If these two classes of objects 
are indeed related, then they should have similar host galaxy 
properties, modulo possible evolutionary effects. In §4.2, we 
found good agreement between the host luminosities of the 1- 
Jy ULIGs and those of infrared-excess quasars (Surace et al. 
2001) and optical quasars from the sample of McLeod & Rieke 
(1994b; McLeod & McLeod 2001), but less so when comparing 
with the sample of radio-quiet and radio-loud quasars of Dun- 
lop et al. (2002). Morphological classification of the hosts in 
ULIGs indicates a preponderance of elliptical-like hosts among 
single-nucleus systems (§4. 1 and §5.2), a result which is at least 
qualitatively consistent with the conclusions of Dunlop et al. 
(2002; note that the conclusions of McLeod & McLeod 2001 
are ambiguous in this respect). We also found that the half- 
light radii of the 1-Jy ULIG hosts are similar to those found 
by McLeod & McLeod (2001), but significantly smaller than 
those measured by Dunlop et al. (2002). The half-light sur- 
face brightnesses of ULIG hosts at R also appear to be slightly 
brighter than those measured by Dunlop et al. (McLeod & 
McLeod 2001 do not tabulate this quantity in their paper). 

It is not clear at present whether the differences found be- 
tween our results and those of Dunlop et al. are significant. 
Many effects may complicate the interpretation of our data. Ex- 
cess /?-band emission from circumnuclear starbursts, known to 
be common in ULIGs (e.g., Veilleux et al. 1995) but less so 
among classical quasars, could explain the brighter fi e in ULIG 
hosts but may have difficulties explaining the ~ 50% difference 
in r e between the ULIG hosts and the quasar hosts of Dunlop 
et al. (2002). As discussed in §5.2, residuals from the AGN 



subtraction may also artificially brighten the half-light surface 
brightnesses and decrease the half-light radii in our sample but 
this effect is not thought to be sufficient to reconcile our re- 
sults with those of Dunlop et al. (2002). The apparent lack 
of agreement between McLeod & McLeod (2001) and Dunlop 
et al. (2002) in terms of the half-light radii suggests that the 
quasar results may depend on sample selection (although these 
two samples have similar mean redshifts, < z > = 0.2, and lu- 
minosities) or the methods of observations (/?-band versus H- 
band HST imaging). On-going adaptive optics imaging studies 
of quasar hosts (e.g., Guyon 2002, in prep.) will increase the 
sample size and will span a broader range in quasar properties; 
the results from these studies should settle this issue within the 
next few years. 

In §4.2, we mentioned the existence in our data of slight 
trends between host luminosities, IRAS 25 ^m/60 ^m colors, 
and nuclear optical spectral types. The most straightforward 
explanation for these trends is that quasar activity generally re- 
quires the mergers of massive (> L* + L*) progenitors. How- 
ever, this scenario has difficulties explaining the fact that the 
trends appear to be stronger at K' than at R. Given the pres- 
ence of circumnuclear starbursts in many of these objects (e.g., 
Veilleux et al. 1995; Kim et al. 1998), it is possible that these 
starbursts contribute significantly to the host luminosities, per- 
haps more so in warm Seyfert galaxies than in the cooler ob- 
jects. If the effect is to be more significant at K' than at R, then 
the starburst age has to be > a few x 10 7 yrs at which point the 
K' luminosity is dominated by young, red supergiants (Bruzual 
& Chariot 1993; Leitherer et al. 1999). Another possible ex- 
planation, also discussed in §4.2, is that the trends are due to 
errors in the PSF subtraction, which are more likely to affect 
the warm highly nucleated Seyferts than the cooler more dif- 
fuse LINERs and H II region galaxies, and are more severe at 
K' than at R. An object-by-object comparison with the //-band 
adaptive-optics results of Surace & Sanders (1999) was carried 
out in §4.2 and seems to indicate a slight luminosity offset for 
the few warm Seyferts in common between the two samples. 
Needless to say, it will be important to expand the sample of 
warm ULIGs observed with adaptive optics or with HST to re- 
fine the measurements of the host luminosities in these objects 
and verify the trends detected in our data. 

6. CONCLUSIONS 

An R- and K'-bmd atlas of the IRAS 1-Jy sample of 118 
ULIGs was presented in a companion paper (Kim, Veilleux, 
& Sanders 2002). The present paper discusses the results from 
the analysis of these images and combines them with the results 
from published spectroscopic studies of ULIGs at optical, near- 
infrared, and mid-infrared wavelengths as well as new Keck 
spectroscopy. The results on the 1-Jy sample are compared with 
those from optical and near-infrared studies of quasars and nor- 
mal ellipticals. The main conclusions are as follows: 

1. All but one object in the 1-Jy sample show signs of a 
strong tidal interaction/merger in the form of distorted or 
double nuclei, tidal tails, bridges, and overlapping disks. 
Interactions involving more than two galaxies are seen 
in only 5 (4%) of the 118 systems. These results confirm 
those of previous studies. 

2. Objects with red R-K' colors and large nuclear-to-total 
luminosity ratios show a tendency to host Seyfert nu- 
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clei, have warm IRAS 25 (im/60 colors and be AGN- 
dominated according to the ISO classification scheme. 

3. Using a classification scheme first proposed by Surace 
(1998) and based on the results of published numerical 
simulations of the mergers of two galaxies, we classi- 
fied the 1-Jy sources according to their overall morphol- 
ogy and apparent nuclear separations. None of the 1-Jy 
sources appears to be in the early stages (first approach 
or first contact) of a merger. Most (56%) of them harbor 
a single disturbed nucleus with and without tidal tails; 
they are therefore in the late stages of a merger. The frac- 
tion of advanced mergers with a single nucleus increases 
above infrared luminosities of 10 12 5 L Q and decreases 
below 10 12 L Q . 

4. The strengths of the H/3 and Mg lb stellar features mea- 
sured in the nuclei of ULIGs are not particularly good 
indicators of the merger phase or epoch of the merger 
event. 

5. All Seyfert Is and most of the Seyfert 2s are advanced 
mergers, either based on their overall morphology or 
their small (< 5 kpc) nuclear separations. A similar re- 
sult is found when we consider the warm objects with 
fio I fm > 0.2 or the AGN-dominated objects based on 
ISO classification. LINERs and H II region-like galaxies 
show no preference between pre-merger and advanced 
merger phases. 

6. The individual galaxies making up the binary systems of 
the 1-Jy sample show a broad distribution in host abso- 
lute magnitudes (luminosities) with a mean of -21.02 ± 
0.76 mag. (0.85 ± g;^ L*) at R and -23.98 ± 1.25 mag. 
(0.90 ± H\ L*) at K', and a luminosity ratio at R or 
K' generally less than ~ 4. The hosts of single-nucleus 
ULIGs have mean absolute magnitudes (luminosities) of 
-21.77 ± 0.92 mag. (1.69 ± g$ L*) at R and -25.03 ± 
0.94 mag. (2.36 ± f |g) at K'. These magnitudes are sim- 
ilar to those found in previous studies of ULIGs, except 
those derived from shallow HST snapshots which un- 
derestimate the contribution from low surface brightness 
features. Correlations are observed between host galaxy 
luminosity and infrared luminosity, IRAS 25 ^m/60 ^m 
color, and optical spectral type. The trend with infrared 
luminosity is due to a redshift bias. There is consider- 
able overlap between the host galaxy luminosity distri- 
bution of single-nucleus ULIGs and that of quasars, al- 
though the hosts of the quasars studied by Dunlop et al. 
(2002) are slightly more luminous on average than the 
1-Jy ULIG hosts. This luminosity shift is not observed 
when comparing with the results of McLeod & McLeod 
(2001) and Surace et al. (2001). The R-K' colors of 
ULIG hosts are similar to those of quasars. 

7. An analysis of the surface brightness profiles of the host 
galaxies in single-nucleus ULIGs reveals that about 35% 
and 2% of the R and K' surface brightness profiles are 
fit adequately by an elliptical-like R l / 4 -law and an ex- 
ponential disk, respectively. Another 38% are equally 
well fit by either an exponential or an elliptical-like pro- 
file. The remainder (26%) of the single-nucleus sources 
cannot be fit with either one of these profiles. Com- 
bining these results, we find that a de Vaucouleurs pro- 



file is an adequate fit to ~ 73% of the single-nucleus 
ULIGs in the 1-Jy sample. These elliptical-like hosts are 
most common in merger remnants with Seyfert 1 nuclei 
(83%), Seyfert 2 characteristics (60%) or mid-infrared 
(ISO) AGN signatures (80%). 

8. The hosts of ULIGs have half-light radii (< r e >=4.80± 
1 .37 kpc at R and < r e > = 3.48 ± 1 .39 kpc at K') which 
are similar to those measured by McLeod & McLeod 
(2001) in quasar hosts, but are significantly smaller than 
the quasar hosts studied by Dunlop et al. (2002). The 
origin of this apparent discrepancy between the two 
quasar datasets is not clear at present. The hosts of 1-Jy 
systems follow with some scatter the [i e - r e relation of 
normal ellipticals (especially if Seyfert Is are excluded 
from the analysis due to possible AGN residuals). The 
distributions of luminosities, /?-band axial ratios, and 
half-light radii in single-nucleus ULIGs are also simi- 
lar to those of normal ellipticals of intermediate lumi- 
nosities. The results at K' are more uncertain because 
of the smaller sample size and possible PSF subtraction 
residuals. Elliptical-like hosts in the 1-Jy sample show a 
broader range of boxiness values than normal ellipticals. 

9. Results #7 and #8 provide strong support to the idea that 
some of the ultraluminous infrared mergers in the 1-Jy 
sample may eventually become intermediate-luminosity 
elliptical galaxies under the condition that they get rid of 
their excess gas or transform this gas into stars. The sig- 
nificant fraction of single-nucleus ULIGs with ambigu- 
ous surface brightness profiles or with large boxiness pa- 
rameters indicate that these objects are still feeling the 
effects of the recent mergers. These results confirm the 
predictions of numerical simulations that violent relax- 
ation is very efficient at the center of the merger but less 
so in the outer regions. 

10. The results from this study are generally consistent with 
the evolutionary scenario in which ULIGs are the re- 
sults of a merger of two gas-rich galaxies which first 
goes through a starburst-dominated pre-merger phase 
when the system is seen as a binary, next reaches a dust- 
enshrouded AGN-dominated merger phase once the two 
nuclei have merged into one, and then finally ends up as 
an optical (post-ULIG) AGN where the host is elliptical- 
like and shows only limited signs of the merger. How- 
ever, our results also indicate that many ULIGs in the 
1-Jy sample may not follow this exact scenario. For in- 
stance, approximately 46% of the 41 advanced mergers 
in the 1-Jy sample show no obvious signs of Seyfert ac- 
tivity, while seven of the 45 pre-mergers already show 
Seyfert 2 (but not Seyfert 1) activity. The possible cor- 
relations between host luminosities, IRAS 25 ^m/60 /zm 
color, and optical spectral types may add another twist 
to the merger scenario, suggesting that merger-induced 
quasar activity may require the merger of massive (> L* 
+ L*) galaxies. These trends may also be due to a bright- 
ening of the circumnuclear starburst in the hosts of the 
warm Seyfert ULIGs. However, one needs to be cau- 
tious when interpreting these apparent trends because of 
the difficulty associated with removing the contribution 
of the central AGN in some of the warm Seyfert ULIGs. 
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Deep near-infrared HST and adaptive-optics imaging of 
ULIGs and quasars will help further clarify the nature of their 
host galaxies. These observing techniques should be applied 
to a large sample of warm ULIGs to refine the PSF subtraction 
in these objects and verify the possible trends with host lumi- 
nosities reported in the present paper. The spatial resolution of 
the data presented here is not sufficient to determine the core 
properties of the elliptical-like ULIG hosts and the bulge-to- 
disk ratio of the galaxies in the pre-merger phase. Both high 
spatial resolution and large dynamic range will be needed to 
fully characterize the host galaxies of these ULIGs. A simi- 
lar program should be carried out on the classical quasars to 
help us understand the apparent inconsistencies between the 
various quasar datasets. Finally, three-dimensional spectro- 
graphs on large telescopes will provide the sensitivity and two- 
dimensional spatial coverage needed to derive detailed kine- 
matic maps in several elliptical-like ULIG hosts. Comparisons 
with data on normal old (low-z) and young (high-z) ellipticals 
will help quantify the similarities and differences between these 
two classes of objects. 
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Appendix A: Notes on Individual Objects 
In this Appendix, we briefly discuss the rationale behind the 
interaction classification of the 1-Jy sample presented in Table 
1 and discussed in §3.4. For this task, we used the R and K' 
images presented in Figure 1 of Paper I and new long-slit spec- 
tra obtained on 2000 March 30-31 and 2001 January 22 - 
23 using Keck II with the Echellette Spectrograph and Imager 
(ESI) in the prism mode. This configuration produces 4000 - 
9000 A spectra with a dispersion that ranges from 50 to 300 km 
s -1 pix" 1 . Accurate wavelength calibration was derived by fit- 
ting a 9th-order Chebychev polynomial to the line positions of 
xenon and mercury-neon lamps observed with the same setup 
as that of the objects. The redshifts derived from these data are 
generally accurate to within ± 0.001 (1 a). This is not accurate 
enough to allow us to make any statement about the relative mo- 
tion of galaxies in a system but it allows us to determine com- 
panionship and group membership. The exposure times were 
typically 10 minutes for each slit position. Systems for which 
we obtained long-slit ESI spectra are indicated by an asterisk. 

• F00091-0738 (Illb): This system is a close pair with ap- 
parent nuclear separation of only 2.1 kpc. The two nuclei are 
more easily visible in the K' image. 

• F00188-0856* (V): This is a single-nucleus system with 
slight distortions in the outer isophotes. The bright source south 
of the nucleus is a star. 

• F00397-1312 (V): This is a single-nucleus system with 
elongated nuclear isophotes and distorted outer isophotes. 

• F00456-2904 (Ilia): The redshift of the north-east compo- 
nent of this pair is unknown. A broad tidal tail is apparent on the 
north-east galaxy and a shorter one extending to the north-west 
is apparent on the main galaxy. 

• F00482-2721 (Illb): This is a relatively close pair with ap- 
parent nuclear separation of 6.7 kpc. The "wisp" to the south- 
west is barely apparent in the K' image and appears to be unre- 
lated to the pair. The redshift of this feature is unknown. 

• F01004-2237 (V): This is a highly nucleated source with 
slight distortions in the outer isophotes (see HST image of 
Surace et al. 1998). 

• F01 166-0844 (Illb): This is a classic looking pair with a 
moderately large nuclear separation of 9.8 kpc and a prominent 
tidal tail extending to the north of the north-west component. 

• F01 199-2307* (Ilia): This is a pair of emission-line galax- 
ies separated by 20.3 kpc. The north-east companion has the 
same redshift as the main component (0.155 versus 0.156, re- 
spectively). A faint tidal tail is visible on the west side of the 
main component, and material also appears to bridge the binary. 

• F01298-0744* (IVb): The IRAS source is the distorted 
object to the north-east. This object harbors two prominent 
tidal tails, signs of a recent merger. The galaxy to the south- 
west shows no sign of interaction and has an absorption-line 
spectrum with a redshift of 0.188 instead of 0.136 for the IRAS 
source. 

• F01355-1814* (Illb): This object is classified as a close 
pair based on the separation (5.8 kpc) between the centroid of 
the bright south-east (SE) component and the fan-like north- 
west (NW) component. The feature to the south-east of this pair 
appears to be a star-forming tidal feature based on the fact that 
it is barely visible at K'. All three features have emission-line 
redshifts of 0.191 to within ± 0.001. Deep multi-band imag- 
ing or spectroscopy will be needed to confirm the nature of the 
faint south-east feature; the existence of an old underlying stel- 
lar population in this object would require changing the classi- 



fication of this system to a triplet. 

• F01494-1845* (IVa): This diffuse single-nucleus source 
shows a broad tidal feature to the north-west. The faint ob- 
ject to the south-east is a star while the elongated source to 
the north is an emission-line galaxy at the same redshift as the 
IRAS source. The compact source to the east of the elongated 
galaxy is a broad-lined quasar with redshift of ~ 1 .62 based on 
the identifications of the C IV A 1549, C III] A 1909, and Mg II 
A2798 features. 

• F01569-2939 (IVa): This highly distorted source harbors 
a narrow tidal tail stretching over 55 kpc. A shorter counter-tail 
is also probably present on the west side of the object. 

• F01572+0009 = Mrk 1014 (IVb): This single-nucleus 
source presents two curved tidal features extending 46 kpc (also 
see Surace et al. 1998). 

• F02021-2103 (IVa): This highly distorted object presents 
two symmetric tidal features. 

• F0241 1+0353* (Illb): This apparently complex system at 
R is a close (7.1 kpc) pair at K'. The north-east (NE) com- 
ponent presents a fan-like tidal feature extending to the east. 
Keck spectroscopy has confirmed the physical interaction be- 
tween the NE and SW components (z = 0.144 and 0.143, re- 
spectively). The compact source near the western edge of the 
/?-band field appears to display a very faint feature extending to 
the south and therefore may not be a star. 

• F02480-3745 (IVa): This single-nucleus object presents 
two symmetric tidal features extending over 27 kpc. 

• F03209-0806 (IVb): This compact object harbors spectac- 
ular tidal tails extending over 61 kpc. 

• F03250+1606* (IVb): This object shows weak tidal fea- 
tures to the north and south-east of the nucleus. The faint com- 
pact object near the eastern edge of the field of view is a star, 
and the same also appears to be true for the other compact ob- 
ject west of the star. The compact object (labelled "G") to the 
south-west has the same redshift (both emission and absorption 
lines) as the IRAS source, but shows no sign of interaction. 

• FZ03521+0028 (Illb): This is a close pair separated by 
only 3.5 kpc. 

• F04074-2801* (IVa): This source is highly distorted. A 
tidal feature is seen extending ~ 28 kpc north of the nu- 
cleus. The two galaxies to the north-west and north-east of the 
IRAS source show no signs of interaction but are at the same 
emission-line redshift as the IRAS source to within the mea- 
surement uncertainties. 

• F04103-2838 (IVb): This compact source presents two 
short curved tidal features extending 8 kpc. 

• F0431 3-1649 (IVa): This is a classic merger remnant with 
two large tidal tails extending 59 kpc. 

• F05020-2941* (IVa): The IRAS source is the distorted 
source with short tidal tails to the south and north-west. The 
bright source to the east is a star. No obvious emission or ab- 
sorption feature was detected in the faint diffuse object between 
the star and the IRAS source. 

• F05024-1941* (IVa): The IRAS source is the distorted 
source near the center of the image. A tidal feature is clearly 
visible on the west side of this object. The east source is a bright 
star. 

• F05156-3024 (IVb): This single-nucleus object presents a 
short tidal feature to the south. 

• F05189-2524 (IVb): This compact galaxy has a faint 
curved tidal feature to the north-west (refer to Surace et al. 1998 
for a high-resolution HST image of this object). 
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• F07599+6508 (IVb): The IRAS source lies to the south- 
east in this field, and presents faint tidal features on the west 
and south sides (more easily visible in the images of Surace 
1998). The other bright source in the field is a star. 

• F08201+2801* (IVa): This object harbors a broad tidal tail 
that extends over 25 kpc. The spectrum of the faint source to the 
north-west presents faint [O II] A3727 emission at the same red- 
shift as the IRAS source, but it is not clear whether the emission 
really originates in the north-west source; it may be contamina- 
tion from the diffuse tidal feature surrounding the IRAS source. 
No obvious emission or absorption feature was detected in the 
southern object. 

• F08474+1813* (V): This apparently complex system com- 
prises the IRAS source to the west and two stars to the south- 
east and north-east of this source. No obvious tidal tail is ap- 
parent in the IRAS source. 

• F08559+1053* (IVb): A diffuse tidal tail seems to emerge 
from the south-west side of the bright galaxy. The two galax- 
ies to the west of the brighter source have the same emission- 
line redshifts as that of the main source to within ± 0.001, but 
show no sign of interaction with the main object. This system 
is therefore classified as a single-nucleus source, and the other 
two galaxies are considered (non-interacting) group members. 

• F08572+3915 (Illb): This is a well-known interacting pair 
separated by 5.7 kpc (e.g., Sanders et al. 1988a,b; Armus et al. 
1990; Surace et al. 1998). 

• F08591+5248* (V): The IRAS source to the east is a com- 
pact object with slightly distorted isophotes but no obvious tidal 
tails. The source to the west is a faint emission-line galaxy at 
the same redshift as the IRAS source (z = 0.157), but it shows 
no sign of interaction. 

• F09039+0503 (IVa): This is a single-nucleus system with 
a prominent curved tidal tail extending 24 kpc to the east. 

• F091 16+0334* (Ilia): The compact source located east of 
the bright Seyfert 2 nucleus is a small absorption-line galaxy at 
the same redshift as the Seyfert galaxy. Tidal features are seen 
to the east and west of the Seyfert galaxy. 

• F09463+8141 (IVa): This source appears to be a pair at R 
but the K' image shows that it is in fact a single nucleus with a 
prominent tidal tail extending out to 76 kpc to the west. 

• F09539+0857 '(V): This compact object has distorted outer 
isophotes but shows no obvious tidal tails. 

• F10035+2740 (IVa): This object appears to be double at R 
but the K' image shows that it is in fact a single nucleus with a 
prominent tail extending out to 54 kpc to the north. The redshift 
of the bright galaxy to the south-west is unknown; it shows no 
obvious sign of interaction with the main object. 

• F10091+4704 (IVa): This highly distorted source presents 
two complex tidal features to the south and north-west of the 
nucleus. The southern feature presents a luminous knot at R 
which is not apparent at K' . 

• F10190+1322 (Illb): This is a close (5.5 kpc) pair with 
weak tidal features. The compact source to the south is star. 

• F10378+1 108 (IVb): This single-nucleus source presents 
broad tidal features extending to the north and south. 

• F10485-1447* (Ilia): Stars complicates the classification 
of this system. The two brightest stars in the field have been 
confirmed spectroscopically. The identification of the third one 
located next to the galaxy labelled "E" is based on its unre- 
solved point spread function. The eastern (absorption-line) and 
western (emission-line) components have similar redshifts (z = 
0.133). Tidal features seem to emanate from the western com- 



ponent and one of them possibly creates a bridge with the east- 
ern component. This interpretation is not unique, however. It is 
possible that the eastern component is not interacting with the 
western component and that the apparent "bridge" is instead a 
second tidal feature from the western component as a result of 
an earlier merger. 

• F10494+4424 (IVb): This object presents two peculiar 
"spikes" to the north which are presumably of tidal origin and 
a fainter diffuse feature to the south-west. 

• F10594+3818 (Illb): This close (4.4 kpc) pair shows 
prominent tidal tails extending out to 90 kpc to the west. The 
bright source to the west is a star. 

• F11028+3130 (IVa): This single-nucleus source shows a 
broad tidal feature to the west. 

• Fl 1095-0238 (IVb): A prominent curved tidal tail is vis- 
ible to the north of this compact source. The high-resolution 
//-band image of Bushouse et al. (2002) reveals an interest- 
ing double core in this object. However, as mentioned in their 
paper, it is not clear whether the double core represents two 
galactic nuclei or a single nucleus bisected by a dust lane. This 
is the reason why this object is classified as a late merger (IVb) 
rather than a close binary (Illb). 

• Fl 1 1 19+3257 {Wb): This compact source shows an inter- 
esting bootlegged tidal feature to the east. 

• Fl 1 1 30-2659 (IVa): This single-nucleus system presents a 
complex tidal feature to the north-west and a shorter feature to 
the south-east. 

• F11180+1623* (Ilia): The IRAS source is the wide (21.8 
kpc) pair of galaxies to the south-east in this field. The emis- 
sion lines in the spectrum of the western component in this pair 
are quite faint. No spectrum was obtained of the two galaxies 
labelled "G". 

• Fl 1223-1244 (Ilia): This source is conservatively classi- 
fied as a very widely separated (88 kpc) pair. The redshift of 
the eastern component is unknown. The signs of interaction in 
the eastern component are very slight and it is possible that the 
tidal features of the western component were caused instead by 
an earlier galaxy merger. 

• Fl 1387+41 16 (V): This compact source shows slight outer 
isophotal distortions. The redshifts of the two faint sources 
north-west of the nucleus are not available. 

• Fl 1506+1331 (IVb): This is a single-nucleus object with 
two symmetric tidal tails to the east and west of the nucleus. 
The nature of the object to the south is unknown, but does not 
appear to participate in an interaction. 

• Fll 582+3020 (V): The IRAS source is the distorted object 
on the east in this field. The compact object to the west is a star, 
and the other two diffuse sources (labelled "G") are galaxies 
with redshifts similar to that of the IRAS source but they show 
no sign of interaction with it. 

• FZ1 1598-01 12 (IVb): The IRAS source is the object in the 
middle of the field with two prominent tidal tails extending out 
to 12 kpc. The redshift of the galaxy to the north is unknown. 

• F12018+1941 (IVb): The IRAS source presents two sym- 
metric tidal tails. 

• Fl 2032+1 707 (Ilia): The K' image of this object reveals 
the presence of two interacting galaxies separated by 12.0 kpc. 

• F 12072-0444 (IVb): This single-nucleus source presents a 
curved tidal tail to the south and a shorter feature to the north- 
west (see high-resolution image of Surace et al. 1998). 

• F121 12+0305 (Illb): This peculiar-looking object is a 
close pair separated by 3.8 kpc. A bright tidal tail is present 
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in the /?-band image. 

• F12127-1412* (Ilia): The members of this widely sepa- 
rated (21.9 kpc) pair have similar emission-line redshifts (z = 
0.133) although the emission lines in the north-east component 
are considerably stronger than those in the south-west compo- 
nent 

• Fl 2265+02 19 = 3C 273 (IVb): This object is highly nu- 
cleated, but shows a thick tidal feature towards the north-east 
and a very faint feature to the north-west (e.g., Tyson, Baum, 
& Kreidl 1982). The faint linear feature emerging to the south- 
west is the well-known optical jet in this system (e.g., Bahcall 
et al. 1995). 

• F12359-0725* (Ilia): This apparent triplet is made of a 
widely separated pair (labeled "N" and "S"; separated by 22.9 
kpc) and of a background galaxy at z = 0.396 to the west. 

• F12447+3721 (IVa): This diffuse object shows broad tidal 
features to the east and north-west. 

• F12540+5708 = Mrk 231 (IVb): This object is made of a 
compact nucleus surrounded by a diffuse tidal complex (e.g., 
Surace et al. 1998). 

• F13106-0922* (IVa): The IRAS source is the highly dis- 
torted source to the north-east in this field with a very long (82 
kpc) tidal tail. The source north of the IRAS source is a star and 
so is the source to the west. 

• F13218+0552* (V): This is a highly nucleated source with 
no obvious tidal tails. The compact object to the south-west 
of the IRAS source appears to be a background absorption-line 
galaxy at z = 0.23 although additional data are needed to con- 
firm this redshift. 

• F 13305-1739* (V): This highly nucleated galaxy shows 
slight isophotal distortions. Our spectrum of the compact 
source to the north-west is not of sufficient quality to determine 
the nature and redshift of this object. 

• F13335-2612 (Illb): The two nuclei in this close (3.2 kpc) 
pair are more easily visible in the K' image. 

• F13342+3932* (IVb): This object is misclassified as an in- 
teracting group by Borne et al. (2000). The IRAS source is in 
the middle of the field. The galaxy to the north-east (labelled 
"G") has an emission-line redshift of 0.180, very similar to that 
of the IRAS source. However, it shows no sign of interaction 
with the IRAS source. The elongated galaxy to the north-west 
lies in the foreground at z = 0.165. No redshift is available for 
the other galaxy to the north-west. 

• Fl 3428+5608 (IVb): This merger remnant harbors a spec- 
tacular tidal tail which extends out 3 1 kpc to the south. 

• F13443+0802* (Tpl) All three objects in the field are 
emission-line galaxies with similar redshifts (z = 0.135). A 
bridge of material may be present between the eastern binary 
and the south-west component. 

• F13451+1232 (Illb): This is a well-known pair separated 
by 4.0 kpc (e.g., Surace et al. 1998; Evans et al. 1999). 

• Fl 3454-2956 (Ilia): A tidal tail appears to emerge south- 
west of the southern component in this widely separated (15.8 
kpc) pair. 

• F13469+5833 (Illb): This close (4.3 kpc) pair shows 
prominent curved tidal tails to the north-east and south-west. 

• F13509+0442* (IVb): The IRAS source is at the center of 
the field and shows a prominent tidal tail which extends 15 kpc 
to the south. The galaxy near the northern edge of the field of 
view is a background emission-line galaxy at z = 0.188. The 
object to the south-west of the IRAS source is a star. 



• F13539+2920 (Illb): This relatively close (7.0 kpc) pair 
presents a strong tidal feature to the north which gives the ap- 
pearance of a third nucleus in the /?-band image (or /-band HST 
image of Cui et al. 2001). This feature is barely apparent in the 
K' image. 

• F 14053-1958 (Illb): This very close (1.5 kpc) pair is 
barely resolved in the K' image and presents two tidal features 
to the north and south. 

• F14060+2919 (IVa): This is a highly distorted source with 
a strong linear tidal feature to the north-east surrounded by a 
diffuse envelope of material. A single nucleus is observed in 
this object at K'. This object is misclassified as multiple (> 2) 
by Cui et al. (2001). 

• F14070+0525 (V): This compact object shows no obvious 
tidal features but has disturbed isophotes. 

• F 14121-0126 (Illb): This moderately close (9.1 kpc) pair 
shows a prominent curved tidal tail to the north-west and a 
fainter feature which seems to emerge from the southern com- 
ponent. The compact source to the south-west appears to be a 
star based on its unresolved point spread function. 

• F14197+0813 (V): A star to the north is superposed on the 
body of this object. There is no evidence for tidal tails in this 
object. 

• F 14202+261 5 (Ilia): The IRAS source is the pair of galax- 
ies to the north. The redshift of the southern galaxy is unknown 
but since it shows no obvious sign of interaction with the pair, 
we classify this system as a pair instead of a triplet (Cui et al. 
2001). 

• F14252-1550 (Illb): This appears to be a moderately close 
(8.3 kpc) pair with a fairly large luminosity ratio of ~ 8 at K'. 

• F14348-1447 (Illb): This is a close (4.8 kpc) binary with 
distinct tidal features to the north and south of the galaxies. The 
compact source near the eastern edge of the field is a star based 
on its unresolved point spread function. 

• F14394+5332 (Tpl) This system is classified as a triplet 
because the eastern component has two close (2.6 kpc) nuclei 
and a prominent tidal feature appears to bridge the eastern and 
western components. However, we have no confirmation that 
the two components lie at the same distance. 

• F 14485-2434 (IVb): The IRAS source lies west of the field 
center. The redshift of the galaxy to the north-east is unknown; 
it shows no sign of tidal disturbance. 

• F15001+1433* (Tpl) This object is classified as a triplet 
because the tidal feature west of the eastern component appears 
to point towards the other two western sources. One of these 
(labelled W in Fig. 1 of Paper I) has the same redshift as the 
main component (z = 0.162). The low signal-to-noise of the 
spectrum for the westernmost (WW) component does not allow 
us to confirm its physical connection with the other galaxies in 
the system. 

• F15043+5754 (Illb): This is a close (5.1 kpc) double with 
short and stubby tidal features. 

• F15130-1958 (IVb): This single-nucleus source has two 
tidal features extending to the north and south over 14 kpc. 

• F15206+3342 (IVb): This is a relatively compact source 
with a distinct tidal feature to the west. HST images of this 
image reveal a complex morphology (Surace et al. 1998). A re- 
cent detailed kinematic study of this object by Arribas & Colina 
(2002) confirms the existence of a single nucleus in this system. 
The source to the north is a star based on its unresolved point 
spread function. 
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• F15225+2350 (IVa): This merger remnant presents a 
prominent tidal tail to the south-west and a fainter one to the 
south-east. 

• F 15327+2340 = Arp 220 (Illb): This system is well known 
to harbor two near-infrared nuclei separated by 0."9 or 0.4 kpc 
(e.g., Scoville et al. 2000), but these two nuclei are not distinct 
at the resolution of our images. 

• Fl 5462-0450* (IVb): The IRAS source presents tidal tails 
extending over 8 kpc. The compact source to the south is a 
spectroscopically confirmed star. 

• F16090-0139 (IVa): This single-nucleus object presents a 
diffuse tidal feature to the north-east. 

• F16156+0146* (Illb): The two objects in this pair have 
emission-line spectra with similar redshifts (z = 0.133). The 
system presents a faint stubby tidal feature to the south and a 
very faint tail to the north-west. A tidal bridge may also be 
connecting the two objects in the pair. The compact source 
near the galaxy pair appears to be a star based on its unresolved 
point spread function, although underlying continuum emission 
from the north-west component makes this statement uncertain. 

• F16300+1558 (V): This is a compact object with distorted 
outer isophotes, but no obvious signs of tidal tails. 

• F 163 3 3+4630 (Ilia): This widely separated (13.1 kpc) bi- 
nary appears to be linked by a tidal bridge. The object to the 
north-east of the pair is a star based on its unresolved point 
spread function. 

• F16468+5200 (Illb): This moderately close (7.8 kpc) pah- 
presents a tidal tail which seems to emerge from the western 
component and extends 14 kpc to the north. 

• F16474+3430 (Illb): This object is similar to 
F14252-1550 in that it is a moderately close (6.5. kpc) pair 
with a fairly large nuclear luminosity ratio of ~ 5.5 at K' (the 
components are too close to each other to determine reliably 
the global luminosity ratio). A prominent tidal tail is visible on 
the west side. 

• F16487+5447 (Illb): This close (5.4 kpc) pair appears to 
be linked by a tidal bridge. 

• Fl 7028+581 7* (Ilia): The western component of this 
widely separated (23.3 kpc) pair presents diffuse tidal features 
and the isophotes of the eastern component are extended in the 
direction of the western companion. Both components have the 
same redshift within the uncertainties (z = 0.106 versus 0.107, 
respectively) 

• Fl 7044+6720 (IVb): The brighest source in the field is a 
star. The IRAS source presents a peculiar tidal feature to the 
north-west of the nucleus. 

• F 1 7068+4027 (Tpl): This source is tentatively classified as 
a triplet because of its resemblance with the other triplet candi- 
date F15001+1433. The main source (labelled "E") in the field 
has a single nucleus but is elongated along the east-west direc- 
tion, i.e. in the same direction as the two fainter western (W 
and WW) galaxies. Note that the redshifts of these two galax- 
ies are unknown and that no obvious signs of interaction are 
visible in these objects. If future observations do not confirm 
the interaction, the IRAS source should be re-classified as a IVa. 

• F17 179+5444* (IVb): This single-nucleus source is sur- 
rounded by what seems to be tidal material extending asymmet- 
rically to the south-west. The compact source to the south-west 
is a star based on its unresolved point spread function and its 
spectrum. 

• F20414-1651* (IVb): The isophotes of this single-nucleus 
object are highly elongated. The spectra of the two compact 



and diffuse sources to the south of the main galaxy are not of 
sufficient quality to allow us determine the nature of these ob- 
jects (a faint Ha feature at the same redshift as the main galaxy 
may be present in the spectrum of the compact object, but this 
feature is very faint and may be due to contaminating emission 
from the main object). These southern sources may be galax- 
ies in interaction with the main source but the evidence is not 
sufficiently convincing to classify this system as an interacting 
pair. 

• F2 1208-05 19 (Ilia): This wide binary is separated by 14.1 
kpc and shows distinct tidal features which extend out to 14 
kpc. 

• F21219-1757 (V): The isophotes of this compact object 
are only slightly distorted. 

• F21329-2346 (IVa): This single-nucleus object presents a 
faint tidal feature to the south and a peculiar bootlegged feature 
to the north which is also presumably of tidal origin. 

• F21477+0502 (Tpl) This is arguably the best case for a 
triple system in the 1-Jy sample. A bridge of material appears 
to connect the western component to the eastern pair. 

• F22088-1831 (Illb): The members of this pair are sepa- 
rated by 4.3 kpc. The compact source to the north-east is a star 
based on its unresolved point spread function. 

• F22206-2715 (Illb): This moderately close (7.7 kpc) bi- 
nary shows two prominent tidal features to the south-east and 
north-west. It is misclassified as a multiple (> 2) system by Cui 
etal. (2001). 

• F22491-1808 (Illb): This close (2.2. kpc) pair shows two 
classic tidal tails to the east and north-west. It is misclassified 
as a multiple (> 2) system by Cui et al. (2001). 

• F22541+0833 (Ilia): This system is tentatively classified a 
widely separated (18.5 kpc) binary. The outer isophotes of the 
south-east galaxy are highly elongated in the direction of the 
north-west companion. Note, however, that the redshift of the 
south-east galaxy is unknown. 

• F23060+0505 (IVb): This single-nucleus system presents 
a diffuse tidal feature to the south-west. 

• F23129+2548 (IVa): The IRAS source lies north in this 
field and present tidal features to the south and west. The red- 
shift of the galaxy to the south is unknown. This galaxy shows 
no sign of interaction with the IRAS object. 

• F23233+281 7 (Isolated): This source shows a spiral-like 
pattern surrounding a very compact nucleus. No obvious tail 
is visible in this object. The point sources directly south and 
north-east of the nucleus are presumed to be bright H II regions 
based on their luminosities and colors. This is the only object 
in the 1-Jy sample with no obvious sign of interaction. 

• F23234+0946 (Illb): This moderately close (7.4 kpc) pair 
shows a hook-like tidal feature to the south-east. The nature of 
the two faint point sources to the south-east and north-west of 
the pair are unknown. Neither one shows signs of interaction 
with the IRAS source. 

• F23327+2913 (Ilia): Each member of this wide (22.7 kpc) 
binary presents tidal features and appears to be linked to each 
other by a bridge of material. 

• F23389+0300 (Illb): This close (4.9 kpc) binary presents 
a short tidal feature to the south. 

• F23498+2423 (Ilia): This wide (12.7 kpc) binary source 
presents prominent tidal tails at R. The compact source west of 
the binary is a star based on its unresolved point spread func- 
tion. 
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FIG. 1 . — Color - magnitude diagrams of the 1-Jy sample, (a) Integrated R-band absolute magnitudes versus R-K'. (b) Integrated K' absolute magnitudes versus 
R-K'. (c) Nuclear (4-kpc diameter) S-band absolute magnitudes versus nuclear R-K'. (d) Nuclear (4-kpc diameter) K' absolute magnitudes versus nuclear R-K'. 
Nuclear quantities are measured within a diameter of 4 kpc. 



FIG. 2. — Trends between the integrated R (left) and K' (right) absolute magnitudes and (a) - (b) the IRAS 25 /im/60 fim colors; (c) - (d) the optical spectral 
types; (e) - (/) the ISO spectral types; (g) - (h) the infrared luminosities. 



FIG. 3. — Trends between the integrated (left) and nuclear (4-kpc; right) R-K' colors and (a) - (b) the optical spectral types; (c) - (d) the IRAS 25 /im/60 fim 
colors. No trends are observed with the infrared luminosities, IRAS 60 fim/ 100 fim colors and the strengths of the H/3 and Mg lb stellar features (not shown here). 



FIG. 4. — Trends between the R (left) and K' (right) compactness values and (a) - (b) the optical spectral types; (c) - (d) the IRAS 25 fim/60 fim colors; (e) - (/) 
the equivalent widths of the stellar Mg lb feature. 



FIG. 5. — Apparent nuclear separations in the 1-Jy sample of galaxies. The distribution is highly peaked at small values but also presents a significant tail at high 
values. The very uncertain separation measured in Fl 1223-1244 (87.9 kpc) is not shown in this figure. 



FIG. 6. — Trends between the apparent nuclear separations and (a) the optical spectral types; (£>) the IRAS 25 /im/60 fim colors. No trends are observed with the 
nuclear values of the equivalent widths of the H/3 and Mg lb stellar absorption features, two traditional indicators of the starburst age (not shown here). 



FIG. 7. — Apparent nuclear separations as a function of infrared luminosities. The great majority of the extreme ULIGs with log [Lfr/Lq] > 10 12 ' 5 Lq are 
single-nucleus systems and therefore have NS < 2.5 kpc. 



FIG. 8. — Total projected tail lengths as a function of (a) the optical spectral types; {b) the infrared luminosities; (c) the IRAS 25 fim/60 fim colors; (d) the IRAS 
60 fimliOO fim colors; (e) the equivalent widths of the H/3 stellar feature; (/) the equivalent widths of the Mg lb stellar feature. No obvious correlations are observed 
with any of these parameters. 



FIG. 9. — Interaction classes for the objects in the 1-Jy sample. See text in §3.4 for a description of each class. The five triplets in the sample are plotted separately 
on the left. The merger of two galaxies is expected to follow the sequence I — » V. Contrary to systems of lower infrared luminosities which often comprise two 
distinct galaxies (classes I, II and III), most (56%) of the objects in 1-Jy sample are single-nucleus systems and are therefore in the later stage of a merger (IV or V). 



FIG. 10. — Positive trends between the interaction class and (a) the infrared luminosity, (b) the optical spectral type, and (c) the IRAS 25 fim/60 fim color. No 
trends are observed with the IRAS 60 fim/ 100 fim colors and the strengths of the H/3 and Mg lb stellar features (not shown here). 



FIG. 11. — Observed surface brightness profiles of (left) R-band images and (right) K' images. These profiles have not been corrected for cosmological dimming. 
The profile of each image is plotted on two different radial scales: a linear scale and ai? 1 / 4 scale. Solid lines in each of the panels are least-square exponential or de 
Vaucouleurs fits to the data. To avoid seeing effects, only the data outside of twice the radius of the seeing disk are considered in the fits. Arrows indicate the size 
of the seeing disk. The value of the reduced chi-square is indicated on the upper right hand corner of each panel. Values much larger than unity indicate poor fits, 
while values much less than unity indicate fits which are poorly constrained. 



FIG. 12. — Fraction of (a) single-nucleus and (b) double-nucleus 1-Jy ULIGs with disk-like (D), elliptical-like (E), elliptical/disk-like (E/D) and ambiguous (A) 
radial profiles. E/D profiles are equally well fitted by an exponential or de Vaucouleurs function, while neither fits galaxies with "ambiguous" profiles. As expected, 
the fraction of ambiguous profiles is considerably larger among binary systems. 



FIG. 13. — Trends between host galaxy profiles and (a) optical spectral types, (b) ISO spectral types, (c) IRAS 25 fim/60 fim colors, and (d) interaction classes. 
Elliptical-like systems show a tendency to host warm AGN and to lie in advanced mergers. The nature of the host does not seem to depend strongly on the infrared 
luminosity, the 60 /im/100 fim ratio, or the strengths of the stellar H/3 and Mg lb features (not shown here). 



FIG. 14. — Absolute R (left) and K' (right) magnitudes of the host galaxies for individual objects in (a) - (b) double-nucleus systems and (c) - (d) single-nucleus 
systems. For comparison, the absolute magnitudes at R and K' of a L* galaxy are -21.2 and -24.1, respectively. 



FIG. 1 5 . — Trends between the S-band absolute magnitude of the host galaxies and (a) the infrared luminosities; (b) the IRAS 25 fim/60 fim colors; (c) the optical 
spectral types. Luminous host galaxies show a tendency to host IR-luminous and warm AGNs. The tendency with infrared luminosity is due to distance-dependent 
effects. 



FIG. 16. — Distribution of the host magnitude differences at (left) R and (right) K' among the binary systems (class III) of the 1-Jy sample. Notice the difference 
in the magnitude scale between the two figures. The magnitude difference is slightly larger at K' than at R, although most systems have a magnitude difference in 
both bands of less than 1.5 or, equivalently, a luminosity ratio of less than 4. 
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FIG. 17. — The nuclear (4-kpc diameter) R — K' colors plotted against the colors of the host galaxies. HII galaxies, LINERs, Seyfert 2s and Seyfert Is are 
represented as stars, open circles, filled circles, and filled squares, respectively. No significant correlation is observed between these two quantities, especially when 
considering only the objects with (R-K') ho!}t < 4. 



FIG. 18. — Seeing-deconvolved half-light radii, r e , of the elliptical-like host galaxies in the 1-Jy sample measured at (a) R and (b) K' . The mean r e of ULIGs at R 
(4.80 ± 1.37 kpc) is significantly smaller than the mean values measured by Dunlop et al. (2002) for the hosts of radio-quiet quasars (7.63 ± 1.11) and radio-loud 
quasars (7.82 ± 0.71 kpc), although there is some overlap. In contrast, the mean r e of ULIGs at K' (3.48 ± 1.39 kpc) is nearly identical to that of the quasars 
measured in the //-band by McLeod & McLeod (2001; 3.39 ± 1.90). The origin of the discrepancy between the two quasar datasets is not known. 



FIG. 19. — Surface brightnesses, fi e , as a function of the half-light radii, r e , measured at (a) R and (b) K' for the elliptical-like host galaxies in the 1-Jy sample. 
The surface brightnesses have been corrected for cosmological dimming and the half-light radii for seeing effects. The hosts of these ULIGs appear to be smaller 
and slightly brighter at R than the hosts of the radio-loud and radio-quiet quasars from the sample of Dunlop et al. (2002), but still follow with some scatter the 
He - r e relation of normal (inactive) ellipticals [represented by the dotted line in panel (a); Hamabe & Kormendy (1987) assuming V — R = 0.5 mag]. At K', the 
interpretation is more difficult because data exist for a fewer number of ULIGs and no quasars. The dash line represents the fit of Hamabe & Kormendy (1987) 
assuming V—K' = 3.7 mag, while the dotted line is the best fit to the A"-band data of Pahre (1999) for normal ellipticals. The hosts of Seyfert 1 ULIGs are noticeably 
brighter (especially at K') than normal galaxies; this may be due to positive residuals from the PSF subtraction. The solid line in Fig. 19i> is the best fit through the 
ULIG data, excluding the Seyfert 1 galaxies. 



FIG. 20. — Radial profiles of the axial ratio in the hosts of single-nucleus ULIGs. The arrow in each panel indicates r = 3 X r e . 



FIG. 21. — Axial ratios measured at 3 X r e in single-nucleus ULIGs. The median and mean of this distribution are 0.77 and 0.75 ± 0.15 (1 <r), respectively. These 
numbers are similar to those measured in normal ellipticals and in quasar hosts. 



FIG. 22. — Radial profiles of the normalized /?-band boxiness parameter, An/a, in the host galaxies of single-nucleus ULIGs. Note the large variations of An/a 
with radius in many of these objects. The arrows in each panel indicate r = 2xr e and 3xr f . 



FIG. 23. — /?-band boxiness values normalized to the semi-major axis at (a) 2xr, and (b) 3xr ( for ULIGs with elliptical-like host galaxies. Both distributions 
are broader than that of normal ellipticals and peak near zero. 
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TABLE 3 
Host Galaxy Magnitudes 
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